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ABSTRACT 


The  gyration  operator  whijh  generates  a  matrix 
from  a  given  matrix  is  defined.  It-  is  shown  that 
the  gyration  is  related  to  the  concept  of  combina¬ 
torial  equivalence.  The  hybrid  matrices  of  net¬ 
work  theory  are  gyration  matrices.  The  ABCD  and 
hybrid  matrices  are  combinatorially  equivalent  to 
their  impedance  and  admittance  matrices. 

■  All  combinatorially  equivalent  matrices  have 
the  same  degree.  The  gyration  operator  also  pre¬ 
serves  the  PP  property  as  well  as  the  rank  of  the 
hermifcian  part. 

The  gyration  operators  form  an  involutary 
Abelian  group.  The  PR  property  is  shown  to  ex¬ 
tend  beyond  the  impedance  and  admittance  matrices 
of  a  passive  network,  and  a  complete  set  of  PR 
matrices  is  given  for  the  description  of  the  net¬ 
work.  The  relationships  between  the  scattering 
and  gyration  operators  are  detailed. 

A  complete  synthesis  procedure,  based  on  the 
gyration  operator  is  given.  Any  PR  immittance 
matrix  can  be  synthesized.  Three  worked  examples’ 


are  included. 


Of  special  significance  is  a  novel  concept 
of  enumerating  the  degree  of  a  rational  matrix 
function. 
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CHAPTER  I 


The  Gyration  Operator 

1.1  INTRODUCTION.  In  this  first  chapter,  we  de¬ 
fine  an  operator  termed  the  gyration.  The  gyra¬ 
tion  operates  on  a  given  matrix  to  produce  another 
matrix  and  is  symbolized  by  the  equation 

T(A)  =  B 

where  A  is  i-  given  matrix,  r  the  gyration  oper¬ 
ator  and  B  i«  the  resultant  matrix,  termed  the 
gyration  of  A. 

This  relationship  is  sufficient  to  make  the 
matrices  A  and  B  combinatorially  equivalent,  a 
term  that  was  coined  by  A.W.  Tucker  p.1]  .  The  im¬ 
pedance,  admittance  chain  and  hybrid  matrices  of 
network  theory  are  all  combinatorially  equivalent. 
The  work  of  A.W.  Tucker  emerged  from  the  linear 
programming  field  and  is  applied  here  to  network 
theory. 

It  is  shown  that  the  gyration  operator  pre¬ 
serves  the  passivity  property  (PR) .  It  is  also 
shown  that  all  combinatorially  equivalent  matrices 
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have  the  same  degree.  (The  degree  measures  the 
complexity  of  a  network. )  A  third  invariant  of 
the  gyration  operator  is  the  rank  of  the  hermitian 
part,  a  property  of  some  significance  in  network 
theory. 

Of  special  inportance  in  this  chapter  is  a 
novel  method  of  enumerating  the  degree  of  a  ration¬ 
al  matrix  function. 

The  relations  between  the  scattering  operator 
and.  the  gyration  operator  are  investigated,  and 
various  theorems  are  detailed. 

1.2  POLE  MULTIPLICITY  OF  A  RATIONAL  MATRIX 

FUNCTION. 

In  a  previous  paper  Duffxn  and  Hazony  (5  ] 
studied  the  properties  of  the  degree  of  a  rational 
matrix  function  F(s).  it  was  brought  out  that  the 
degree  may  be  defined  in  several  equivalent  ways. 
One  of  these  ways  concerned  the  poles  of  the  minors 
of  F  (s)  . 

Definition.  Let  F ( s )  be  a  matrix  whose  ele¬ 
ments  are  rational  functions  of  the  complex  vari¬ 
able  s.  Let  k  be  the  number  of  distinct  poles 


that  occur  in  the  matrix  elements  (the  pole  at 
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inf inity  is  counted) .  Then  the  degree  may  be  de- 
fined  to  be 


(1) 


6F(s)  =  h1  +  h2  +  ...  +  h^. 


Here  h .  is  the  maximum  multiplicity  with 
th 

which  the  j  pole  appears  in  the  minors  of  F(s). 

By  "minor"  is  meant  the  determinant  of  F  or 
the  determinant  of  any  square  submatrix  of  F. 

(in  addition  if  the  submatrix  has  no  rows  or  col¬ 
umns  we  define  the  empty  minor  to  have  the  value 
1.  This  convention  will  be  employed  in  a  lemma 
to  follow. )  F  need  not  necessarily  be  square. 

It  is  convenient  to  have  a  notation  for  the 
multiplicity  of  a  pole  at  a  given  point.  Thus.,  if 
£  denotes  a  complex  number  o::  the  point  at  in¬ 
finity  let  6^.F(s)  be  defined  as  the  maximum  multi¬ 
plicity  of  the  pole  at  s  =  (  of  any  minor  of  the 
matrix  F(s).  We  term  6^F  the  multiplicity  of  F(s) 
at  With  this  notation  we  can  describe  Relation 
1  so  as  to  define  the  degree  of  F  in  the  form 


(2) 


6F(s)  =  Z  6 *F (s) 

c  ' 


Here  the  summation  is  over  all  points  of  the  com- 
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plex  plani*  including  the  point  at  infinity.  When 
the  operator  6^,  is  applied  to  a  scalar  function, 
it  will  not  conflict  with  the  usual  meaning  of 
multiplicity. 

To  derive  properties  of  the  multiplicity 
operator  6^,  the  following  lemma  is  very  useful. 
Lemma  1 .  Let  M  and  N  be  n  by  n  matrices.  Then 
the  determinant  of  their  sum  is 


(3) 


M  +  N|  = 


Z 

i  =  0 


Ml  JN 


Here  the  sum  is  over  all  minors  | M | ^  of  M  multi¬ 
plied  by  the  algebraic  complementary  minor  [ N | 1 
of  N.  Also  |m|q  =  1/  l N| °  =  | N[ ,  |M|k  »  |M|, 
and  |n[  =1,  where  k  is  the  total  number  of  minor 
determinants  in  an  n  n  matrix.  The  proof  of 
this  result  is  given  in  a  previous  paper  [5]. 

Lemma  2.  If  6^G  =  0,  then  6^ (F  +  G)  =  6^F. 

Proof.  Let  H  =  F  +  G?  and  suppose  H',  F',  and  G' 
denote  square  submatrices  of  H,  F,  and  G.  Let 
6^(F'|  -  6^(F);  but  6^  |F"  | <6^ (F)  if  the  submatrix 
F"  has  fewer  rows  than  F' .  Lemma  1  is  now  applied 
with  M  =  F'  ar.d  H  =  G'  .  Then  it  follows 


) 


f* 


that  one  term  on  the  right  side  of  Relation  3  is 
| F  *  I  alone  while  all  other  terms  contain  minors 
of  F  with  fewer  rows  than  F' .  Thus  the  term  |F'| 
has  a  greater  multiplicity  than  any  other  term. 
Hence, 

6^(F'  +  G')  =  6^|F'  +  G' |  =  6^|F'[  =  6^F 

This  shows  that  d^H  -  d^F.  But  F  =  H  -  G,  so  by 
a  symmetrical  argument  it  follows  thi  6^F  -  d^H. 
Thus  6^F  =  6^.H,  and  the  proof  is  complete. 

Lemma  3  to  follow  is  stated  here  for  complete¬ 
ness;  it  will  be  used  in  a  later  paper. 

Lemma  3.  Let  4>  (s)  be  a  scalar  function  and  F(s) 
be  an  n  by  n  matrix  function.  Then 

d^  ( 4>F )  ^  a^F  +  na^a  . 

Proof.  Let  H  =  $F,  and  suppose  that  d^H  =  d^jH’l 
for  some  m  by  m  minor  H' .  Then 

a^H  =  6^| 4>p •  |  =  a^ <4>mj f*  | )  *  a^m  +  a^|F*  [  , 

but  a^4>m  =  raa^0  -  na <t>  and  a^|F'|  -  d^F,  so  the 
proof  is  complete. 
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1.3  GYRATION  OF  A  MATRIX. 

Given  a  square  matrix  A  with  matrix  elements 
a^j.  If  a^  0,  we  define  a  matrix  B  with  matrix 
elements 


(4) 


bi4  =  "ai5/an 

V  “  aiii/,an 

V  =  VrauiaiC/an 


£  =  2,3,... 
|i  ~  2,3,... 


The  matrix  B  is  termed  the  gyration  of  A  about 
pivot  a^.  This  operation  of  forming  a  gyration 
may  be  denoted  by  r.  Thus, 


B  =  F  (A)  . 


It  is  clear  then  from  the  above  definition  that 
A  =  r  (B) .  In  other  words  FF  =  I.  The  term 
gyration  was  suggested  by  the  fact,  that  a  gyrator 
transforms  an  impedance  matrix  A  into  B. 

Let  xl  be  an  arbitrary  column  vector.  Then 
A  defines  a  transformation  yj =  Axl  .  In  a  three 
dimensional  case  this  stands  for 


* 


T 
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so  we  can  eliminate  x^  from  the  right  side  of 
these  equations  and  obtain 


Note  that  the  last  matrix  is  precisely  B,  the  gyra¬ 
tion  of  A.  It  follows  that  f  be  defined  as  a  trans¬ 


position  of  x^  and  y^  and  symbolized  as 
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(5)  r(xx  x2  x3  yi  y2  y3)  =  (yx  x2  x3  x3  y2  y3)  . 

A  similar  definition  applies  to  a  gyration  of  A 

about  any  other  diagonal  element  a^  as  pivot. 

This  gyration  may  be  denoted  as  I\.  Then  it  is 

x 

seen  from  the  transposition  property  of  that 
if  i,  j,  k  be  a  permutation  of  the  integers 
1,  2,  3/ 

a*1  =  ri(rj[rk(A)i) 

This  gives  a  rapid  way  (closely  related  to 
Gaussian  elimination)  of  computing  the  inverse 
of  a  matrix.  For  example,  if 


A 
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and  the  last  matrix  is  A  \  Then  it  may  be  said 
that  the  gyration  of  a  matrix  is  a  partial  in¬ 
verse  of  a  matrix. 

It  will  be  recognized  that  the  gyration 
operation  is  precisely  that  of  forming  the  hybrid 
matrix,  well  known  in  network  theory  from  the 
analysis  of  series  parallel  conf igurationsjVj. 

A  transformation  more  general  than  gyration 
has  been  introduced  by  A.W.  Tucker  [11].  Tucker 
calls  his  relation  between  matrices  combinatorial 
equivalence.  it  follows  from  his  Theorem  7  that 
two  matrices  are  combinatorially  equivalent  if 
and  only  if  it  is  possible  to  pass  from  one  to 
the  other  by  a  finite  succession  of  elementary 
operations  of  the  following  three  types: 

1.  interchanging  of  any  two  rows 

2.  interchanging  of  any  two  columns 

3.  A  gyration 

It  is  clear  from  this  characterization  and 
the  transposition  property  of  gyrations  that  a 
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comb  ina  tor  ial  equivalence  can  be  described  as  a 
permutation  of  the  sequence  (x^  x^  y^  y ^  ^^3) 
such  as  (yi  x3  x2  y2  x^  . 

A  basic  property  of  combinatorial  equi¬ 
valence  is  given  in  Tucker’s  Theorem  3  which 
follows • 

Let  {a}  denote  the  set  of  square  submatrices 
of  A  of  all  orders  (including  an  empty  submatrix 
<t>  of  order  zero)  .  Let  { (3}  denote  the  like  set  for 
B  where  B  is  combinatorially  equivalent  to  A.  Then 
there  is  a  one-to-one  correspondence  (3-«— ►a  between 
{(3}  ana  such  that  corresponding  subdeterminants 
| Pj  and  |a[  are  proportional  within  sign.  Speci¬ 
fically 

Oi  =  ±M/|a*|  , 

where  a  =  a*  corresponds  to  (3  =  <j>  (taking  |  <t> |  =  1)  . 
The  nonsingular  matrix  a*  is  called  the  pivot  of 
the  transformation  of  A  into  B. 

Tucker  says  that  his  work  on  combinatorial 
equivalence  was  suggested  by  the  well  known  simplex 
method  devised  by  G.  B.  Dantzig  for  the  solution 
of  linear  program  problems,  we  wish  to  show  in 
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this  chapter  that  combinatorial  equivalence  has 
application  in  the  entirely  different  field  of 
network  synthesis.  It  is  worth  noting  that  ideas 
similar  to  combinatorial  equivalence  have  been 
applied  to  network  algebra  problems  by  Bott  and 
Duff  in  [2]  . 

We  now  apply  our  results  on  the  degree  of 
rational  matrix  functions  together  with  Tucker's 
Theorem  3. 


Theorem  1.  Let  A  be  a  rational  matrix  function  of 
the  complex  variable  s  and  suppose  B  is  combina¬ 


torially  equivalent  to  A.  Then  B  is  a  rational 


matrix  function  of  the  same  degree  as  A. 


Proof.  Let  p  denote  |a*|,  the  pivot.  Then  accord¬ 
ing  to  Tucker's  theorem  | £|  =  |a|p-1.  Let  x  be 

-1 

neither  a  pole  of  p  or  p  ,  and  apply  the  multi¬ 


plicity  operator  6  .  It  follows  that 

X 

6x1  Pi  =  6x(|a[p"1)  =6x|a[.  Then 


6  A  =  6  B 
x  x 


Let  y  be  a  pole  of  p.  Then  it  follows  that 
6  |  pj  -  6  |  aj  -  6  P,  and  so  by  the  definition  of 

y  y  t 


1 
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the  multiplicity  operator 


(5B-6A-6P. 

y  y  y 

Let  z  be  a  pole  of  p_1.  If  6  A  >  0,  let  ot  be  the 

z 

submatrix  so  that  6z|a  |  =  6ZA.  But  if  6ZA  =  0,  let 
a  be  the  empty  matrix,  so  jet  [  =1.  Then  in  either 
case 

6zlpi  =  6zla|  +  6ZP  1  =  6ZA  +  6ZP  1. 


It.  follows  that 


6zB 


6zA  +  V 


-1 


The  degree  of  B  is  defined  as  6B  =  2  6  B  , 

b 

S 

where  the  summation  is  over  all  finite  values  and 
also  the  point  at  infinity,  clearly  the  points 
x,  y,  and  z  are  distinct  but  together  include  all 
s  points.  Thus 


6B=  2  6  B  +  X  6  B  +  2  6  B  . 

X  X  z  Z  y  y 


Substituting  the  relations  just  found  for  6  B, 
6yB,  and  6zB  gives 

6B  —  2  6  A  +  2  6  A  +  X  5  A  -  X  6  P  +  Z  6ZP 

x  y  *  z  Y  J  z 


Ik 

& 
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The  last  two  sums  cancel  since  6p  =  6p  \  Thus 
it  follows  that  6B  -  6A.  But  combinatorial 
equivalence  is  a  symmetric  relationship  so 
6A  -  6B.  This  proves  that  6A  =  6B. 

1.4  CHARACTERIZATION  OF  GYRATION  MATRICES. 

Let  A  be  an  n  by  n  matrix  of  real  numbers, 
let  x|  be  a  vector,  and  let  yj  =  Axf.  The  scalar 
product  of  the  vectors  xj  and  yj  is  denoted  by 
3T  yj  and  defined  as 

n 

*  Yl  =  2  x.y.  • 

1  1  1 

If  xAxJ  is  positive  for  every  real  vector  /  IjT, 
then  A  is  said  to  be  positive  definite.  This 
differs  from  the  standard  definition  in  that  we 
are  not  requiring  A  to  be  symmetric.  We  propose 
to  extend  this  definition  in  a  natural  way  when 
A  has  complex  matrix  elements.  Let  "x  be  an  arbi¬ 
trary  complex  vector  and  let  ~x*  be  the  complex 
conjugate.  Then  if  x*Axf  is  in  the  right  half¬ 
plane  or  on  the  j-axis,  we  shall  say  that  A  is 
is  right  definite.  This  condition  may  be  stated 
as  an  inequality 


Re  "x*Axt  -  0  . 
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Note  that  this  condition  is  less  restrictive  than 
the  condition  x^AxJ  0,  which  is  one  of  the  de¬ 
fining  conditions  for  non-negative  definite  her- 
mitian  matrices. 

Theorem  2.  The  gyration  of  a  right  definite  matrix 
A  is  a  right  definite  matrix  B. 

Proof.  it  is  assumed,  of  course,  that  the  pivot, 
say  does  not  vanish.  Let  uj  be  an  arbitrary 

vector  and  let  v|  =  BuJ .  Let  the  vector  x|  be 
defined  as  follows; 

X1  =  V1,X2  =  U2'  X3  =  U3'  ***'  Xn  =  Un  * 

If  y|  =  AxJ  ,  then  by  the  property  of  a  gyration 

Y1  =  ul'  *2  "  v2'  ^  =  V  Thus 

u*v|  =  x*y(  -  x^*y^  +  x.^* 

Taking  the  real  part  gives 

Re  u*v|  =  Re  x*y{  or 
Re  u*Bu»  =  Re  x*Ax{  -  0  . 


This  completes  the  proof. 
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Definition.  A  matrix  function  F(s)  of  the  complex 
variable  s  is  said  to  be  positive  real  (PR)  if  I 
and  II  are  satisfied: 

I.  The  matrix  elements  f  ^  (s)  are  rational 
functions  of  s  with  real  coefficients. 

II.  For  any  choice  of  complex  numbers 

xl"**'  xn 
n  n 

Re  2  2  f.  .x.x.*  Ss  0  for  Re  s  ^  0  . 

1  i  13  3  i  - 

Condition  II  states  that  F  is  a  right  definite 
matrix  for  s  in  the  right  half  plane.  Of  course 
Condition  II  is  not  meaningful  at  a  pole. 

Let  f(s)  -  x*Fx{  .  Then,  as  is  shown  in 
Appendix  1,  Condition  II  is  equivalent  to  the  con¬ 
dition  : 

II  a.  Re  f(s)  ^0  for  Re  s  =  0 

II  b.  f(s)  has  no  poles  for  Re  s>  0. 

II  c.  For  Re  s  =  0  poles  of  f(s)  are  simple 
and  have  non-negative  residues. 

The  impedance  (or  admittance)  matrix  of  a 
passive  network  is  shown  Appendix  1  to  be  positive 
real.  if  the  network  has  no  gyrators  it  is  re¬ 
ciprocal;  i.e.,  the  matrix  F  satisfies  the 
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symmetry  conditions: 

III.  f..  *  f ..  io£  3  =  x'  2>  •••' 

1  J  J  1 

we  are  concerned  in  this  chapter  mainly  with 
reciprocal  networks.  Thus,  let  G  be  the  gyration 
of  a  matrix  F,  which  satisfies  I,  II,  and  III  and 
f  ^  does  not  vanish  identically.  Then  G  may  be 
regarded  as  the  hybrid  matrix  of  a  reciprocal 
network.  G  satisfies  I.  By  Theorem  2,  G  satis¬ 
fies  II,  and  by  Equation  4  which  defines  a 
gyration,  the  symmetry  Condition  III  is  replaced 
by: 


=  2,  3, 


Condition  III1  may  be  termed  (first)  hybrid 


Note  that  if  G  has  hybrid  symmetry  and  g-^fs) 
does  not  vanish  identically,  then  r (G)  has  regular 
symmetry.  Thus  we  may  state  the  following: 


Theorem  3.  conditions  I,  II,  and  ill  are  neces- 


sary  and  sufficient  that  G  be  the 


H 


tion  of  the  impedance  matrix  of  a  passive  net- 


rovidea  gi;L(s)  does  not 


work  without 
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vanish  identically. 

Proof,  This  is  a  consequence  of  Theorem  2  and  the 
above  definitions. 

A  positive  real  matrix  F  will  be  said  to  be 
IPR  if 


n  n 

IV.  Re  2  2f..x.x.*-0  for  Re  s  =  0 

1  1  13  3  1 

As  is  well  known,  an  L  C  network  (iue.,  a  net¬ 
work  without  resistors  or  gyrators)  is  character¬ 
ized  by  having  an  impedance  (or  admittance)  matrix 
F  satisfying  I,  II,  III,  IV. 

Theorem  4.  Conditions  I,  II,  III1,  and  IV  are 
necessary  and  sufficient  that  G  be  the  (first) 
gyration  of  the  impedance  matrix  of  an  L  C  net¬ 
work,  provided  g^^s)  does  not  vanish  identically. 
Proof.  Let  F  =  r (G)  and  let  Re  s  =  0.  Then  given 
a  vector  IT  we  have  shown  in  the  proof  of  Theorem  2, 
that  there  is  a  vector  ~x  so  that  Re  x*Fxl  =  Re  u*Gul 
Conversely,  given  an  arbitrary  IT  there  exists  a  IT, 
which  satisfies  this  equation.  Then  Re  u*Gul  =  0 
for  all  ~u  if  and  only  if  Re  x*Fx(  =  0  for  all  3T. 

But  the  latter  condition  characterizes  L  C  networks, 


y 
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so  the  proof  is  complete. 

It  is  known  that  matrix  elements  of  an  IPR 
impedance  matrix  are  odd  functions  of  s.  Likewise 
it  can  be  seen  from  Equation  4  that  if  G  is  a 


hybrid  IPR  matrix,  g^  and  are  odd  functions 
while  g  ,  and  g.„  are  even  functions. 

V-  14 

If  a  matrix  G  satisfies  Conditions  I,  II, 
and  III'1',  then  G  is  a  positive  real  matrix 


function  with  hybrid  symmetry.  However,  if  g^(s) 
vanishes  identically,  then  it  is  not  the  gyration 


about  the  first  pivot  of  a  symmetric  positive  real 


matrix. 


1.5  REMOVAL  OF  IMAGINARY  AXIS  POLES. 

Let  G(s)  be  the  hybrid  matrix  of  a  passive 
reciprocal  network,  and  suppose  that  G  has  poles 
on  the  j  axis.  According  to  condition  II  c,  the 
function 


n  n 

g(s)  =  Z  Z  g..(s)x.x  * 
1  1  13  1  3 


can  have  only  simple  poles  on  the  j  axis.  Since 

the  numbers  x^  are  arbitrary,  it  is  seen  that 

g.  .  (s)  can  have  at  most  simple  poles  on  the  j  axis. 


V 
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The  proof  goes  as  follows: 

First  take  x.  =  6-u  where  6..  is  the 

1  lh  lh 

Kronecker  delta.  Then  we  see  that  g(s)  =  g^(s) 
and  so  g^fs)  can  have  only  simple  poles  (by  II  c) . 
Next,  take  x^  =  6^  +  an<^  so 

g(S)  =  gu(s)  +  ghh(s)  -  j[glh(s)  -  ghl(s)] 

Making  use  of  the  hybrid  symmetry  gives 

g(s)  =  gil(s)  +  ghIl(s)  +  2jghi(s) 


Since  g,  g^,  and  g^h  have  at  most  first-order 

poles  on  the  j  axis,  it  follows  that  g^  has  a 

first  order  pole  at  most. 

Next  take  x.  =  +  6-u  for  k  ^  2,  h  ^  2. 

i  lk  in 

Then  similar  considerations  show  that  g^  has  a 
first-order  pole  at  most  and  the  proof  is  complete. 

Thus  consider  a  pole  at  the  point  s  =  sq  =  jco0 
on  the  j-axis  and  suppose  sq  ^  0.  Then 


a . 


g  .  .  ( s )  =  - ^ —  +  g '  .  .  ( s ) 

yi;j  '  '  s  -  so  v  13 ' 


where  a. .  is  the  residue  constant  and  g'  .  .  (s)  is 
13  il 


T 
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bounded  at  s  =  s  .  Let  s  -  s  =  £  be  a  number  in 

o  o 

the  right  half  plane.  Then  for  £  of  small  absolute 
value,  it  is  seen  the  first  term  on  the  right  is 
dominant.  consequently  by  Condition  III'*’  it  is 
seen  that  the  matrix  A  has  hybrid  symmetry. 

Then  for  arbitrary  x 


x*Gxi 


x*Axl 

-  + 

s  -  s 

o 


x*G "  xl 


The  first  term  on  the  right  is  dominant 

for  |  £ |  small.  Hence,  it  follows  from  Condition 

II  c  that  "x*Ax|  -  0.  in  particular,  x*Ax|  is  real. 

By  Condition  I  it  follows  that  there  is  a 

pole  at  s  =  s  *  =  -jm  and 
o  o 


(S) 


where  g"  .  .  (s)  is  bounded  at  s  =  +  jo)  .  Let 
j  ij  o 


Of  course,  A*  satisfies  the  same  condition  as  A. 
Thus  P  satisfies  Conditions  I  and  III*’.  However, 
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since  x*Ax|  and  x*A*x|  are  both  real  non-negative 
it  follows  that  the  quadratic  form  >c*Pxl  is  a  sum 
of  poles  on  the  j  axis  with  real  non-negative 
residues.  Thus  P  satisfies  both  II  and  IV,  and  so 
P  is  hybrid  I PR.  The  matrix  G"  =  G  -  P,  so 
clearly  I  and  ill'*'  are  satisfied.  II  a  is  satisfied 
because  Re  x*G"x|  =  Re  x*Gx|  on  the  j  axis.  II  b 
is  satisfied  because  neither  G  nor  P  has  poles  in 
the  right  half  plane.  Note  that  G"  has  no  poles 
at  s  =  +  ja>  so  at  any  other  j-axis  point,  residue 
residue  (x*G"x| )  =  residue (k*Gx| ) .  Hence,  G"  satis¬ 
fies  II  c.  This  t. roves  that  G"  is  a  hybrid  PR 
matrix. 

Theorem  5.  If  G  is  a  hybrid  PR  matrix  then 

G  =  P  +  Q 

where  P  is  a  hybrid  I  PR  matrix  and  Q  is  a  hybrid 
PR  matrix  without  poles  on  the  j-axis. 

Proof.  if  the  only  poles  of  G  are  at  +  ja)Q,  then 
the  above  argument  gives  the  proof  with  Q  =  G" . 
Otherwise  the  process  is  repeated  on  G",  etc.  It 
is  seen  that  a  pole  at  zero  or  infinity  can  be  re¬ 
moved  by  a  similar  method. 


r 


1.6  SCATTERING  MATRICES  AND  HYBRID  MATRICES f 3] 

Let  A  be  a  matrix  such  that  (A  +  I)  ^  exists. 
Then  the  scattering  operator  S  is  defined  as 

S  (A)  =  (A  -  I)  (A  +  I)"1  . 

Theorem  6.  A  is  a  right  definite  matrix  if  and 
only  if  T  =  I  -  S S  is  a  non-negative  definite 
hermitian  matrix  (where  S^*  is  the  conjugate 
transpose  of  S,  and  S  stands  for  S (A) ) . 

Proof.  Given  x|  let  y{  =  (A  +  I)  ^xl  so 
x|  =  (A  +  I)yj  .  Then 

x*Tx|  =  x*xj  -  x*St*Sx| 

=  y*(At*+I) (A+I) yj  -  y*(At*-l) (A-I)yj 
=  2y*Ay|  +  2yA*y*{ 

Then  x*TxJ  =  4  Re  (y*Ay|  )  and  the  proof  is  complete. 

It  is  desired  to  find  the  scattering  matrix 
S (B)  of  the  gyration  matrix  B  -  T (A) .  It  is 
sufficient  to  consider  the  case  of  three-by-three 
matrices.  First  it  is  clear  that 
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B  - 


This  can  be  written  in  the  form  (omitting  zeros) 
(6)  B  = 


Moreover  B  +  I  = 


and  B  -  I  = 


Since  (A  +  I)  is  assumed  to  exist,  it  follows 
that  (B  +  I)  ^  exists.  This  shows  that 
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S  (B)  =  (B  -  I)  (B  +  I) 


(A  -  I)  (A  +  I)  .  This  proves 


the  following: 

Theorem  7.  Let  A  be  a  matrix  such  that  ^  0 

and  (A  +  I)  ^  exists.  Let  B  be  the  gyration  of  A. 

Then  S(B)  =  JS (A) ,  where  S (A)  is  the  scattering 

matrix  of  A,  S (B)  is  the  scattering  matrix  of  B, 

and  J  is  the  identity  except  =  -  1. 

2 

Now  xt  is  noted  that  J  =  I  so  we  have 
Corollary:  T(B)  =  T(A)  so  B  is  PR  if  and  only  if 
A  PR. 

Now  suppose  that  A  is  a  rational  matrix 
function. 

Theorem  8.  If  A  is  a  rational  matrix  function  and 
(A  +  I)  ^  exists,  then  the  degree  of  A  and  the 
degree  of  the  scattering  matrix  S (A)  are  equal. 
Pr°of.  Since  S (A)  =  I  -  2 (A  +  I)"1,  we  see  that 
6S  =  6  (A  +  I)-1  =  6  (A  +  I)  =  6A. 

Corollary;  The  degree  of  A  is  equal  to  the 
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proof.  S(B)  =  JS (A)  and  since  J  is  a  nonsingular 
constant  matrix,  6S(B)  =  6S (A) .  Thus  6B  =  6A. 
Theorem  9.  Let  p,  q  be  any  two  complex  numoers 
for  which  B (p) ,  B(q)  are  defined,  and  where 
B  =  T (A) .  Then 

B(P)  +  B(q)T  =  M(q)T[A(p)  +  A(q)T]M(p), 

/  1  -a.l2<P>  -al3(P>\ 


where  M(p)  = 


7^(py  all(p)  an(p) 


Proof.  By  Equation  6  above, 


B  =  CD  ^  where 


Hence 


B(P)  +  B(q)T  =  C  (p)  D_1  (p)  +  [C(q)D_1(q)]  T 
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=  D"1(q)T  [D(q)TC(p)  +  C(q)TD(p)]  d"1  (p) 

Now 

D(q)TC(p)  +  C(q)TD(p)  =  A(p)  +  A(q)T 

also 
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also  to  the  paraherrait ian  part  of  the  matrix  A; 
i.e,,  -i[A(s)  +  A(-s)T]f  subject  to  a  similar  set 
of  assumptions  as  in  the  above  corollary. 


CHAPTER  II 


The  Gyration  Operator  in  Networks. 

2.1  INTRODUCTION.  In  this  chapter,  the  r  operator 
is  applied  to  network  synthesis  and  a  complete  pro¬ 
cedure  is  given  for  the  realization  of  symmetric 

PR  impedance  (or  admittance)  matrices. 

2.2  HYBRID  MATRICES  IN  NETWORK  ANALYSIS. 

Suopose  we  are  given  a  network  whose  imped¬ 
ance  matrix  is  Z(s).  Then  the  currents  and  volt¬ 
ages  are  related  by 

Vi  =  Z  l| 

Let  H  =  r (Z)  where  it  is  assumed  for  the  moment 
that  r (Z)  is  defined.  Partitioning  vj  and  il 
as  follows 
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we  have  by  Equation  5, 


rf  V1  V2^  il  I2^  =  f  il  V2^  V1  I2^  J 


li  \ 

/v,\ 

1 

=  H 

— 

i‘»7 

Consider  the  impedance  Z'  and  Z"  connected  as 

shown  in  Figure  1  with  ports  1  in  parallel  and 

* 

ports  2  through  n  in  series. 

By  definition 

V'J  =  Z'  I'l 
V"t  =  Z"  I"i 

Letting  H'  =  r(Z')  and  H"  =  r(Z")  we  have 


It  is  assumed  that  this  connection  does  not 
change  the  immittance  properties  of  either  net¬ 
work.  This  can  always  be  ensured  by  the  appropri¬ 
ate  use  of  isolation  transformers. 
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lv] 

X 

II 

: 

< 

i  ~  _ 

1  CM 

1  > 

iv‘/ 

1  1 

/  Vl« ' 

— 

X 

II 

— 

I  1 

1  V/ 

But  v1*  =  v1"  and  l2'j  =  i2''{  . 

Hence, 


/v  +  V'l 

=  (H'  4-  HM) 

!  < 
i 

■ 

\V>  +  V1] 

i  v 

We  have  thus  proved  the  following: 

Theorem  10.  Let  two  networks  be  connected  with 
ports  1  in  parallel  and  ports  2  through  n  in 
series.  Then  the  hybrid  matrix  for  the  over-all 
network  equals  the  sum  of  the  hybrid  matrices 
for  the  individual  networks. 

2.3  SYNTHESIS  PROCEDURE. 

The  network  synthesis  problem  consists  of 
associating  a  passive  electrical  network  with  a 
prescribed  PR  matrix  function.  The  basic  approach 


* 


T 
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is  to  extract  simple  sections  from  the  matrix 
that  can  be  realized  by  inspection,  thereby  re¬ 
ducing  the  degree  of  the  matrix.  When  the  degree 
is  reduced  to  zero,  the  process  terminates.  These 
simple  sections  are  extracted  in  such  a  way  as  to 
ensure  that  what  remains  is  still  PR  and  hence 
realizable. 

Assume  that  we  are  given  a  symmetric  PR 
impedance  matrix  Z  to  be  synthesized.  As  a  first 
step,  all  obvious  imaginary  axis  poles  are  re¬ 
moved.  These  are  readily  synthesized.  Once  this 
is  completed,  various  schemes  are  used  to  induce 
further  imaginary  axis  poles,  which  are  removed 
and  synthesized.  One  of  the  methods  of  inducing 
these  further  poles  in  the  scalar  case  is  due  to 
Brune  [12]. 

In  what  follows,  we  shall  use  the  r  operator 
to  give  a  new  extension  of  the  Brune  synthesis  to 
n  port.  Such  extensions  already  have  been  made 
by  a  number  of  people.  •  Brockway  McMillan  [8]  pro¬ 
ceeded  with  the  removal  of  a  certain  amount  of 
resistance  from  each  port  of  Z  until  the  even- 
part  matrix  Ev  Z(s)  =  -^[Z(s)  +  Z(-s)]  becomes 
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singular  at  some  point  on  the  j  axis,  say  ju^. 
B.D.H.  Tellegen  [10]  has  shown  in  his  method  that 

A 

Ev  Z  can  be  made  singular  by  the  removal  of  re¬ 
sistance  from  just  one  port. 

Once  Ev  Z(s)  is  singular,  the  odd  part  then 
is  appropriately  modified  at  jo)Q  resulting  in 
Z(jo)0)  being  singular.  Hence,  Z  ( jooo)  ~ has  a 
pole  that  can  be  removed.  V.  Belevitch  [1]  con¬ 
siders  the  Brune  methods  of  the  above  two  and 

» 

extends  their  results  to  nonreciprocal  networks. 
R.W.  Newcomb  [9]  considers  the  nonreciprocal  case 
with  additional  detail.  (The  methods  of  Gewertz 
and  Oono  [6,  p.  276],  both  non-Brune,  consist  of 
removal  of  j-axis  poles  accompanied  by  successive 
matrix  inversions.  A  similar  method  given  by 
Duff in  [4]  shows  that  network  synthesis  can  be 
viewed  as  a  purely  algebraic  process). 

Our  extension  of  the  Brune  method  differs 
from  the  above  in  that  it  is  not  necessary  at  any 
stage  to  invert  a  matrix.  As  will  be  shown,  the 
method  is  minimal  in  the  sense  of  the  following 
theorem,  stated  by  Tellegen  [10]: 


Z5*  PT'"’ 

* 


— i** 


Minimal  Theorem.  A  positive  real  symmetric  matrix 


function  of  degree  d  can  be  synthesized  as  the 
impedance  matrix  of  a  network  having  a  total  of  d 
reactive  elements.  By  a  reactive  element  we  mean 
either  an  inductor  or  a  capacitor.  An  ideal  trans¬ 
former  is  not  regarded  as  a  reactive  element. 

One  of  the  results  of  the  work  to  follow  is 
a  new  proof  of  the  Minimal  Theorem. 

Following  the  above  procedures,  we  assume 
that  Ev  Z(ja>0)  is  singular.  it  is  shown  in  [6] 
that  there  exists  a  real  constant  matrix  A  such 
that 

M  A 

Z  =  AT  Z  A 


has  Re  z11(j(uQ)  =  0. 

(in  order  to  compensate  for  the  application 
of  this  congruence  transformation,  we  also  apply 
the  inverse  congruence  transformation,  which  may 
be  realized  by  ideal  transformers  as  shown  in  [6]). 

If  z^^  is  identically  zero  then  it  is  shown 
in  [6]  that  the  entire  first  row  and  first  column 
of  Z  is  identically  zero.  Synthesis  then  resumes 
on  the  rest  of  Z,  We  can  thus  assume  without  loss 
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of  generality  that  z^  is  not  identically  zero. 

Suppose  Re  z^  is  identically  zero  on  the 
j-axis.  Then  z^  is  IPR,  and  since  it  cannot 
have  j-axis  poles  (these  have  all  been  removed) 
it  must  be  that  z^  is  identically  zero.  Since 
this  has  already  been  ruled  out,  we  may  take  it 
that  Re  is  not  identically  zero  on  the  j-axis. 

Finally,  suppose  Z  is  identically  singular. 
Then  it  is  shown  in  [6]  that  there  exists  a  con- 
gruence  transformation  such  that  CT  Z  C  has  its 
entire  first  row  and  first  column  identically  zero. 
As  before,  we  may  assume  without  loss  of  generality 
that  this  is  not  the  case,  i.e.,  Z  is  not  identi¬ 
cally  singular. 

Assume  then  that  Z  is  an  n  x  n  PR  impedance 
matrix  without  j-axis  poles,  not  identically 
singular,  z ^  is  not  identically  zero,  Re  z^  is 
not  identically  zero  on  the  j-axis  and  Re  z^  =  0 
at  jo>0. 

Following  the  classical  Brune  tradition  [12] 
we  now  add  a  scalar  b  =  sL  or  1/sc  to  z^  so  that 
Zn  +  b  is  zero  at  ja)Q.  If  cuQ  =  0  or  oo  then 
Re  ( jai0)  =  0  implies  that  z^ ^  ( jcoQ )  =  0  and  no 


rk  rpsp 

4*  " 

i*' 


7" 
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impedance  b  is  required. 


Then  z1i(j^D)  =  0  but  z^  is  not  identically 
zero.  Since  Z  has  no  j-axis  poles,  it  follows 
that  Z  has  only  the  j-axis  pole  possibly  due  to 
b.  Hence  if  0  u)Q  <  oo  (  then  Z  has  at  most  one 
j-axis  pole  in  only  z^  and  this  pole  is  either 
at  0  or  oo  ,  depending  on  b.  if  =  0  or  oo ,  Z 
has  no  j-axis  poles. 

We  are  now  in  a  position  to  reduce  the  degree 
of  Z  by  the  removal  of  a  lossless  section. 

Theorem  11.  Let  Z(s)  be  an  n  x  n  PR  impedance 
matrix  a)  not  identically  singular,  b)  with¬ 
out  j-axis  poles,  c)  with  Re  z ^  not  identically 
zero  on  the  j-axis  and  d)  Re  z^  =  0  aj:  jto  . 

\  01  i  0i 

scalar  b  is  so  chosen  that  =  0* 

Then  Z  may  be  decomposed  in  the  series  para¬ 


llel  manner  shown  in  Figure  1,  where 
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1.  Z'  is  PR  and  6Z '  =  6Z  -  6Z" . 

2.  Z"  jls  ipr  and  ,pf  degree  1  ^>r  2,  and 

3.  Z"  contains  at  most  1  inductor,  1  capac¬ 
itor  and  1  n-port  ideal  transformer  (1 
core,  n  windings) . 

(These  are  referred  to  as  propositions  1,2,  and  3 
respectively) .  proof  is  deferred. 

This  theorem  is  the  basis  of  the  procedure.  By 
repeated  application  of  this  threorem  together 
with  Theorem  12  to  follow  we  eventually  reduce 

A 

the  degree  of  Z  to  zero  at  .;hich  point  the  iter¬ 
ation  terminates. 

To  apply  this  theorem  at  each  cycle,  we  may 
have  to  add  a  reactance  to  z-j.1*  At  t^ie  c°mPle~ 

tion  of  the  cycle  we  therefore  have  to  subtract 
this  reactance.  It  turns  out,  just  as  in  the 
scalar  case,  that  this  subsequent  negative  ele¬ 
ment  can  be  incorporated  in  a  perfectly  coupled 
transformer. 

Theorem  12.  Let  Z,  Z  and  b  be  as  defined  in 
Theorem  11.  Then  after  the  application  of  that 
theorem  to  split  Z,  it  is  always  possible  to  in¬ 
corporate  the  compensating  negative  reactor  into 
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a  perfectly  coupled  transformer. 

The  proof  is  again  deferred.  The  entire 
synthesis  procedure  is  displayed  in  flow  chart 
form  in  Figure  2. 

Example  1 

For  the  sake  of  clarity,  we  choose  an  example 
that  is  already  tailored  appropriately  for  direct 
application  of  Theorem  11.  Let  Z(s)  be  given  by 


Then  r (Z)  - 


s2+s+l 

T  " 

s  +  1 


\ 


S2+  1 


2(s+l) 

“2  , 
s  +  s  +  1 
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+  1 


\ 


s  +  1 


S  +  S  +  1 


! 


-  p  +  o 


r(p) 


=s  2't 


r(Q)  = 


°  \ 
s  +  ?tt j 


z* 


Hence,  Z  is  realized  by  the  network  of  Figure  3. 
Example  2.  Consider  a  problem  given  in  [6,  p.  275]. 


Z(s) 


2s  +  4s  +  5 

7 

2s  +  s  +  1 


■ 


s  +  2  s  +  2 

7 

2s  +  s  +  1 


s2  +  2s  +  2 
__ 

2s  +  s  +  1 


2 

s  +  s  +  1 
— - 

2s  +  s  +  1 


/ 


-jy- 


Ev  Z 


and 


Thus 


Thus, 


/4  2 

4s  +  8s  + 


yZs4  +  3s2  +  2 


2s  + 


2s4  + 


4  2 

4s  +  3s^  +  .1 


det  Ev  Z  =  (2s2  +  1)  • 

4  2^ 

(4s  +  3s^  +  1) 

det  Ev  Z  is  already  zero  at  +j,  +  j '/2//2 . 

Choosing  0)  =  1  we  have 

o 

h  A 

Re  Z(j)  = 

v  V 

using 


v 


4 


I 
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we  have  Re (A  Z(j)A)11  =  0.  Hence 


s  +  s  +  2 


(s  +  2s  +  3) 


z  -  a_za  -  V“(s  +  28  +  3) 

T  “  ^ 


2s  +  4s  +  5 


2s  +  s  +  1 


and  at  s  =  j,  z^  =  -j, 
Thus 


Z  =  Z  4- 


2(s2  4-  1)  (s  +  1)  ~(s2  4-  2s  4-  3) 


_  \-(s~  4-  2s  4-  3) 


2s  4-  4s  4-  5 


2s  4*  S  4*  1 


has  a  j-axis  zero  in  z^. 

Forming  f  (Z) ,  we  obtain 


-41- 


Since  z ^  was  zero  at  s  =  j,  r (Z)  has  a  pole  at 

s  =  j. 

The  pole  matrix  is 


which  is  of  degree  2.  P  can  be  removed  from  F(Z)  to 
give 
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zm  =  r  (p)  = 


Note  that  the  degree  of  Z"  is  also  2. 

Removing  P  from  r (Z)  leaves  Q  =  r (Z)  -  P. 

Then 


/2s  +  2  -1 

z*  =  r  (q)  =  | 

.  2 

4s  +  s  +  2 
__ 

4s  +  2  s  +  2 

Since  we  added  s  to  port  1  of  Z,  we  must  now 
subtract  it.  Also,  to  get  Z  we  require  the  in¬ 
verse  congruence  transformer  at  the  input.  Note 
that  Z'  and  Z"  posses  series  inductors  in  ports 
1,  which  can  be  combined  with  -s  as  shown  in 
Figure  4. 

A 

The  final  realization  of  Z(s)  is  given  in 
Figure  5. 

Note:  The  Brune  transformer  could  be  regarded  as 

an  ideal  transformer  and  one  inductor.  Thus  four 


\ 

/ 
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reactive  elements  are  required,  which  is  in  ac¬ 
cord  with  the  fact  that  6Z  =  4. 

The  dots  on  the  transformers  have  the  follow¬ 
ing  significance,  when  progressing  along  the 
winding  from  the  dotted  terminal,  each  winding 
encircles  the  core  in  the  same  sense.  Thus,  if 
dots  are  at  the  same  ends  of  the  windings,  the 
coupling  term  in  the  matrix  is  positive;  if  at 
opposite  ends,  negative. 

Proof  of  Theorems  11  and  12.  Before  commencing 
with  the  proof  of  Theorem  11,  it  is  necessary  to 
make  some  comments  concerning  notation,  while 
the  operator  r  was  completely  defined  by  Equations 
4,  it  will  be  convenient  in  what  foil ^ws  to  em¬ 
ploy  the  following  approach. 

Let  Z  be  partitioned  into 


where 
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is  a  scalar, 

2f^0  is  a  row-vector  of  n  -  1  elements, 

is  a  column-vector  of  n  -  1  elements, 

and 


Z22  isann-lxn-1  matrix. 

Regarding  two  n-vectors  as  n  x  1  matrices, 
we  have 


Xiy2 


Xiy3 


Xny2 


x  y_ 
nJ  3 


x  y 
nJn 


which  is  a  rank  1,  n  x  n  matrix.  Letting  0”  be 
a  null  row-vector  and  0  a  null  matrix  it  can  be 
verified  that 


r(z)  =  (zu)-1 

M 

(1  i 

-Z^  + 

1°  i 

1 

1 

!  * 

1 

1  Z2l'i 

^  0} 

1 

i  Z22  / 

is  consistent  with  the  definition  of  F(Z)  given  in 
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Equation  4. 

Proof  of  Theorem  11. 

Proposition  1. 

Form  T  (Z) .  Since  z^  has  discrete  zeros  at 

+  jco  ,  T  (Z)  has  poles  there.  By  Theorem  2  or 

the  corollary  to  Theorem  7,  r (Z)  is  hybrid  PR. 

Hence  by  Theorem  5,  f  (IS)  can  be  decomposed  into 

a  hybrid  IPR  matrix  p,  formed  by  the  removal  of 

the  poles  at  +  ja>o  from  f(Z)  and  a  matrix  Q  « 

r (Z)  -  P,  which  is  hybrid  PR  without  j-axis  pole 

It  is  shown  in  the  proof  of  Proposition  3 

« 

below  that  p^ ^  is  not  identically  zero.  Hence 
Z"  =  T (P)  exists. 

Suppose  is  identically  zero.  Since  p^ 

is  IPR,  it  follows  that  Refp^  +  q.^)  =  0  every¬ 
where  on  the  j-axis.  Hence  by  the  corollary  to 
Theorem  9,  z^  =  r  (P-j^  +  q1?)  is  also  IPR.  But 
this  contradicts  assumption  c)  of  the  statement 
of  this  theorem.  Hence  q11  is  not  identically 
zero  and  so  Z‘  =  r (Q)  exists. 

P  and  Q  share  no  common  poles.  Hence  it 
follows  by  Equation  1  that  6 (P  +  Q)  =  6P  +  6Q. 
(For  a  detailed  proof,  see  Theorem  7  of  [5]). 
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But  by  Theorem  1,  the  degree  is  invariant  under 
a  gyration,  so 

$Z  =  <5P  +  6Q  -  6Z'  +  6Z"  . 


Assume  first  that  co0  ^  o  or  oo .  Then  the  hybrid 
IPR  pole  matrix  removable  from  f  (Z)  is  given  by 
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Writing  V 


/vn  i 

7i2\ 

i 

_ t  _ 

_ 

i 

\V2l'  i 

1 

<  1 

to  1 

we  have  by  Equation  9, 


s  -  3m 


s  =  jm 


say,  which  is  real  and  positive  since  is  PR. 


V12  “  -(s  '  ja,o>  T. 


B  »  3m 


=  -jba  by  Equation  7 


=  ~j0  say. 


Z2ll 

V21|  -  (s  -  jm0>  — 


|S  =  3o>, 


jba|  by  Equation  7 


=  . 


Finally 


-<•  -  n>>  — ^ 


<Z21>  )Z*12 


s  =  JO), 


=  ba(  a-. 
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Hence 


(10)  P 


/b 


2s 


i  2 

s  +  cdo 


4>  i 


Z"  can  be  realised  by  Figure  6. 

When  (0o  =  0  or  oo,  the  above  analysis  leads 
to  Z"  being  a  matrix  of  degree  1  and  consisting 
solely  of  a  shunt  reactor  across  port  1,  with 
short  circuits  at  the  remaining  ports. 

This  completes  the  proof  of  Theorem  11. 
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Proof  of  Theorem  12. 

First  we  require  the  following: 

Lemma  4.  Let  a)Q  /  0  o£  oo.  Then  the  matrix  Z" 
of  Theorem  11  can  be  realized  either  by  a  series 
inductor  and  an  ideal  capacitive  transformer  (as 
proved  previously)  or  by  a  series  capacitor  and  an 
ideal  inductive  transformer. 

Proof:  Observe  that  for  sT  a  real  constant  vector, 

/*  :  ? 

both  G  «=[ - 1 - 

and  -G  are  hybrid  PR  by  Definition  I,  II,  and  III1. 
It  is  worth  noting  that  they  are  realizable  with 
one  transformer  and  one  gyrator. 

We  recall  that  in  the  proof  of  Theorem  11, 
the  matrix  r (Z)  was  split  into  P  and  Q.  The 
situation  is  not  changed  if  we  realize  instead 
P  +  G  and  Q  -  G.  Both  are  hybrid  PR  and  hence 
realizable.  Moreover,  6 (P  +  G)  =  6P  and  6 (Q  -  G) 
6Q  by  Lemma  2, 
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Now,  if  we  let  Q  *=  2(3/0^,  then  using  P  defined  in 
Equation  10, 


P  +  G  = 


s  -f  o> 


2s23  I 


(S  +0)  )(Jd 

'  o  '  o 


2  — 


(s2  +  Oi  2)0i 

'  o  '  o 


(s  +  a>o*)b 


r  (p  +  g)  « 


2  2 
S  +  CD 

_ _o_ 

2sb 


0)  b 
o 


£_  £ 
to  b 
o 


2s  £j_£ 

b 


which  is  realized  in  Figure  7.  This  proves  the 
lemma. 

We  proceed  with  the  proof  of  Theorem  12. 

(A  scalar  version  is  found  in  [12]). 


Consider  first  the  case  where  for  s  *=  jo>  , 


then 


211(*do)  “  -^o 


(y>0)  . 
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Then  We  form  Z  with  elements 


zij ■  ?ij +  ^s6ii6jr 

where  6  means  the  Kronecker  delta. 

But  =  2^  +  7s  has  a  zero  at  +  jcoQ  and  a  pole 
8  -  00.  The  pole  matrix  P  has  a  one-one  element 

2  2  *v 

2sb/(s  +  o>o  ) ,  which  was  removed  from  l/(z ^  +  7s) , 
leaving  as  the  one-one  element  of  Q  the  PR  function: 


f(s) 


s2  +  rn  ^  -  27bs2  -  2bsz 


11 


(s2  +  u>q2)  (zxl  +  ys) 


Recall  that  z^fs)  i-3  PR  and  has  no  pole  at  00. 

The  inverse  of  f(s)  is  the  one-one  element 
of  Z'  which  yields  a  series  pole  of  value 
37/(1  -  27b).  Note  that  7/(1  -  27b)  is  real  and 
positive  since  it  is  the  residue  at  a  j-axis  pole 
removed  from  a  PR  function.  The  first  realiza¬ 
tion  of  this  cycle  is  given  in  Figure  8. 

The  element  -7s  compensates  for  the  7s 
added  to  z.^.  The  above  realization  can  be 
given  the  equivalent  network  shown  in  Figure  9. 

In  case  2^(ju>o)  *=  j  (7/0^)  we  add  7/s. 

Then  taking  account  of  Lemma  4  above,  the  final 
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realization  of  the  Brune  cycle  becomes  that  of 
Figure  10. 

The  upper  part  of  the  diagram  shows  a  "capaci¬ 
tive  transformer".  A  two-port  capacitive  trans¬ 
former  is  equivalent  to  a  capacitor  and  an  ideal 
transformer  (see  [6,  p.  114]). 

This  completes  the  proof  of  Theorem  12. 


wmmzmmmmmmm 


mmmamsmmmmm. 
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CHAPTER  III 


Generalizing  the  Gyration 

3.1  INTRODUCTION.  In  this  chapter  the  definition 
and  properties  of  the  gyration  are  extended.  It 
is  shown  that  sets  of  gyrations  can  be  specified 
which  form  Abelian  groups  of  order  2n.  The  general¬ 
ized  gyration  provides  an  extended  basis  for  the 
description  of  network  parameters. 


3.2  EXTENSION  OP  THE  T  OPERATOR 

Suppose  Yl  «  A  xl  where  A  is  an  n  x  n  matrix, 
xl  and  Y\  are  n-vectors. 

Partitioning  as  follows 


M 

1  i 

/A  n 

i 

!  A12  ^ 
1 

jxA 

!  “cNl 
!  * 

a 

^A21 

I 

i 

1 

S  A22  J 

w 

where  A  is  n  x  n,  A^  is  r  x  r,  and 

r-vectors,  we  find  by  solving  for  Xjl  that 
-1 

A^^  exists: 


are 

if 
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This  interchange  of  and  Y ^  is  precisely 
the  r  operator  when  applied  simultaneously  to  the 
first  r  pivots  of  A.  We  say  then  that 


) 


then 


MTYl  «  M  MTXj 
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Let  MtY(  =  Y'l  ,  =  X'l  ,  and  MTA  M  —  A*  . 

Then  Y'»  =  A'X'I  will  be  a  rearrangement  of  the 
equation  Yl  =  AXl  ,  in  which  any  r  pivots  can 
be  made  to  be  the  first  r  by  an  appropriate  choice 
of  M.  The  r-fold  gyration  on  the  selected  set 
can  now  be  performed  according  to  the  definition 
given  in  Equation  (11)  which  is  thus  completely 
general. 

Without  loss  of  generality  let  r  =  2,  n  =  3. 
Then  from  the  standpoint  of  the  vectors  X|  and  Yj.  , 
Equation  .‘(11)  is  equivalent  to: 

ri,2*Xl  X2  x3  yi  y2  y3*  =  (yl  y2  x3  X1  X2  y3> 

and  it  is  thus  clear  that 

(12)  r1/2(A)  =r1[r2(A)]  =r2[r1(A)j  . 

Of  course  for  r  =  1,  Equation  (11)  is  consistent 
with  the  definition  of  given  in  Equations  (4). 

It  follows  from  Equation  (12)  that  Theorems 
1  and  2  are  true  for  the  r-fold  gyration  opera¬ 
tion. 

Suppose  F  is  a  PR  matrix  and  let  G  be  an 
r-fold  gyration  of  F.  Then  Theorem  3  generalizes 


■SS5,:  .„-se  x:i3t3S!!fflsSBSi'?$,S6RJJE33SSS;5;55: 


„  "7  _ 


to  the  following: 

Theorem  13.  Conditions  I  and  II  are  necessary 
and  sufficient  that  G  be  the  r-fold  gyration  of 
the  impedance  matrix  of  a  passive  network,  pro¬ 
vided  that  in  the  parcitioning  of  G,  the  r  x  r 
submatrix  G^  is  not  identically  singular. 

Theorem  4  can  likewise  be  generalized  as  follows. 
Theorem  14.  conditions i/  II  and  IV  are  necessary 
and  sufficient  that  G  be  the  r-fold  gyration  of 
the  impedance  matrix  of  a  network  without  resistors 
provided  that  in  the  partitioning  of  G,  G^1  is  not 
identically  singular. 

It  is  important  to  note  that  in  the  above 
generalizations  and  in  the  generalizations  to 
follow  all  symmetry  conditions  have  been  completely 
relaxed.  The  theorems  will  thus  be  fully  appli¬ 
cable  to  PR  matrices  in  general,  and  will  not  be 
restricted  to  the  symmetric  case. 


3.3  REMOVAL  OF  IMAGINARY  AXIS  POLES  FROM  NON¬ 
RECIPROCAL  MATRICES 

Suppose  F  is  a  PR  matrix,  not  necessarily 


symmetric.  If  f^  or  det 


ft 


\f21  f22 ) 


f^\is  zero  at 


some  discrete  point  on  the  j-axis  then  it  is  clear 
that  the  gyration  matrices  r^(F)  or  2^  will- 
have  poles  on  the  j-axis.  Let  these  gyration 
matrices  be  called  G.  Then  by  reasoning  as  in 
section  1.5,  we  have  the  following  generalization. 
Theorem  15.  Let  G  be  the  £-fold  gyration  of  a 
PR  matrix.  Then 


G  =  P  +  Q 

where  P  _is  IPR  and  Q  _is  PR  without  poles  on  the 
j-axis. 

3,4  GROUP  PROPERTIES 

For  the  case  n  =  3,  let 


Subject  to  an  appropriate  gyration,  each  of  the 

pairs  (yr  x^ ,  (y^  x2) ,  (y^  x3)  can  be  made  to 

lie  in  the  state  (y^  x^)  as  in  the  above  equation 

or  in  the  gyrated  state  (x^  y^) .  Thus  there  are 
3 

precisely  2  or  8  distinct  gyrations  which  can  be 
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defined  for  a  3  x  3  matrix.  In  case  of  an  n  x  n 
matrix  there  are  2n.  .yhe  existence  of  all  of  the 
2n  gyration  matrices  depends  on  the  principal  sub¬ 
matrices  of  the  matrix  .  A.  The  existence  of  a 
given  yyration  matrix  follows  only  when  a  corres[- 
ponding  principal  submatrix  of  A  is  non-singular. 
Definition.  The  operation  of  following  one 
gyration  by  a  sedond  is  called  a  cascade. 

rm£rn(A)J  9enera^es  a  cascaded  pair  of  operators 

r  r. 

m  n 

Theorem  16.  Let  A  be  an  n  x  n  matrix  for  which 
all  of  the  2n  possible  distinct  gyration  matrices 
exist.  Then  the  gyration  operators  form  an  in- 
volutary  Abelian  group  of  order  2n  over  the  cas¬ 
cade  operation. 

Proof. 


ri.[r  j  (a)  3  =  (a) 

and  so  the  set  is  closed  over  the  cascade  opera¬ 
tion. 

Let  i  j  k  be  integers  which  lie  between  0  and 
n.  Define  rQ  as  the  gyration  which  interchanges 
of  t!hfs  (y^  x^) .  Then 


none 


and  so  the  set  is  associative  over  the  cascade 
operation. 

Let  m  b^  any  sequence  of  integers  between 
0  and  n,  without  repeating. 

Then  ro[rm(A)]  =  rm(A)  and  so  rQ  is  the 
identity  element  of  the  set. 

Since  r  fr  (A) 1  =  A  it  follows  that  every 
element  is  its  own  inverse  and  so  the  set  forms  a 
group  with  the  involutary  property. 

Finally  rm[Tn (A) ]  «  rn[Tm(A)]  and  so  the 
group  is  Abelian. 

This  completes  the  proof. 

The  group  can  be  symbolically  displayed  by 
the  use  ofr  a*  hypercube  of  dimension  n. 

Let  n  =  3;  then  the  cube  is  3-dimensional 
and  can  be  thought  of  as  lying,  in  the  first 
hyper-octant  of  an  orthogonal  triplet  of  axes 
a,b,c.  (See  Figure  34  ), 

The  gyration  matrices  can  be  associated  with 
the  corners  of  the  cube  as  follows. 

Suppose  m  is  an  ordered  string  of  integers 
without  repeats  chosen  from  1,2,3  or  n  =  0. 
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Then  r^CA)  is  a  gyration  matrix  which  is  one 
of  possibly  8  distinct  matrices.  If  m  includes  i, 
then  the  i-th  pair  (y^  x^)  is  interchanged  and  so 
the  i-th  coordinate  is  given  the  value  1  in  the 
a,b,c  space.  The  following  table  shov/s  the  re¬ 
lationship  between  rm(A)  and  its  coordinates  in 
the  a,b,c  system. 


TABLE  I 


GYRATION  MATRIX 

COORDINATES 

a  b  c 

r0(A> 

0  0  0 

r^A) 

10  0 

r2(A) 

0  10 

r3(A) 

0  0  1 

rl,2<A> 

110 

ri,3(A) 

10  1 

r2,  3  (A) 

0  11 

rl/2,3(A) 

111 

Figure  '.34  shows  the  location  of  rm(A) 


and  its  coordinates. 
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The  operators  themselves  can  be  associated 
with  spatial  directions  in  the  a,b,c  space. 

r1  will  change  any  matrix  from  its  given 
form  to  a  new  form  along  an  edge  parallel  to  the 
a-axis. 

r2  operates  parallel  to  the  b-axis  and 
parallel  to  the  c-axis. 

As  an  example  r2(A)  has  position  coordinates 
(0,0,1).  The  rL  operation  gives 

ri[r  3  <*)  1  =  r1;3(A) 

which  has  coordinates  (1,0,1)  and  is  given  by  a 
u\ove  along  the  "a"  direction. 

Crossing  the  diagonal  of  a  cube  face  is 
associated  with  a  2-fold  gyration.  There  are 
three  2-fold  gyrations  associated  with  the  three 
planar  senses  of  the  cube  faces,  crossing  a  major 
diagonal  of  the  cube  is  associated  with  the  3-fold 
gyration.  The  3-fold  gyration  is  equivalent  to 
matrix  inversion  in  this  case,  and  so  each  of  the 
gyration  matrices  is  faced  by  its  inverse  at  the 
opposite  end  of  its  corresponding  major  diagonal. 

Starting  with  a  given  matrix  A,  we  may  or  may 
not  be  able  to  perform  a'l  of  the  2n  gyrations 
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possible,  depending  on  the  existence  of  ftIr  Thus 
the  group  may  or  may  not  be  complete. 

Theorem  17.  For  a  given  matrix  A  the  group  of 
operators  is  complete  if  and  only  if  every  princi¬ 
pal  submatrix  of  A  is  non-singular. 

Proof.  Suppose  the  group  of  operators  is  complete. 
Then  it  is  possible  starting  from  A,  .to  move  either 
along  an  edge,  across  the  diagonal  of  a  face  or 
across  a  major  diagonal  of  the  cubo  and  thereby  to 
arrive  at  any  other  corner  of  the  cube. 

Hence,  every  principal  submatrix  of  A  must 
be  non- singular. 

Conversly  suppose  every  principal  submatrix 
of  A  is  non- singular.  Then  by  Equation  11  every 
one  of  the  2n  operators  is  defined. 

This  completes  the  proof. 

3.5  APPLICATION  TO  NETWORK  THEORY 

Suppose  a  passive  network  gives  rise  to  an 
impedance  matrix  Z(s).  Then  Z  is  PR.  Let  it  be 
of  degree  d  and  suppose  that  at  s  =  sq  the  rank  of 
its  hermitian  part  is  r.  If  all  of  the  principal 
minors  of  Z(so)  are  non-zero,  thpn  from  Theorem 
17,  every  cjyration  matrix  is  defined  at  sQ.  It 
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follows  that  a  complete  hypercube  can  be  drawn 
whose  corners  are  associated  with  Z(s  its  in- 
verse  Y(sq)  and  all  of  the  gyration  matrices 
rjzlsjl  .  Each  of  these  matrices  is  PR,  of 
degree  d  and  has  a  hermitian  part  of  rank  m,  a 
result  which  follows  from  Theorems  1, 2  and  the 
corollary  to  Theorem  9. 

If  any  of  the  principal  submatrices  are 
singular  at  sq,  but  not  identically  so,  then  the 
appropriate  gyration  matrix  has  a  pole  at  sQ.  (It 
is  this  property  which  underlies  all  of  the  syn¬ 
thesis  theory  developed  in  Chapters  2  and  4) . 

If  any  of  the  principal  subraatrices  are 
identically  singular,  then  that  gyration  matrix  is 
not  defined  and  the  cube  is  incomplete. 

Consider  the  following  examples. 

Example  1 

2 

A  scalar  z  =  s  t  1  . 

s 

The  hypercube  is  1-dimensional  and  has  two 
corners,  one  at  the  orgin  and  one  at  the  point  1, 
together  with  the  joining  edge.  It  is  complete. 
(See  figure  35) 
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At  the  origin  is  z  with  poles  at  s  =  0  and 
s  =  oo,  zeroes  at  +j  and  a  hermitian  part  which 
is  zero  on  the  j-axis.  z  is  PR  and  of  degree  2. 

At  the  point  1  lies  y  =  z”^  which  has  poles 
at  j-j,  zeros  at  0  and  oo.  y  has  a  zero  hermitian 
part  on  the  j-axis,  is  PR  and  of  degree  2. 

Example  2. 

A  i 

Z  =  s 

l1  1 

Z  is  IPR,  degree  1,  and  the  2x2  matrix  Z 
is  identically  singular  but  its  1x1  principal  sub¬ 
matrices  are  non-singular,  although  zero  at  s  =  0, 
The  2-dimensional  hypercube  is  thus  degenerate 

and  is  shown  in  Figure  36.  it  contains  only  two 

edges  and  the  three  points  (0,0),  (1,0)  and  (0,1). 

At  (0,0)  lies  Z.  At  (1,0)  lies 

/l  .  \ 

r1(z)  = 
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which  has  a  pole  at  s  =  0,  is  IPR  and  of  degree  1. 
At  (0,1)  lies 

r  ‘i 

r2(z)  - 


with  the  same  properties  as  Z. 

Z”1  «  2^  * 8  not  defined. 

Example  3. 

A  unity  turns  ratio,  ideal  transformer  is  certainly 
PR  but  has  no  impedance  or  admittance  matrix.  It 
is  defined  instead  by  the  hybrid  equation 


The  inverse  of  the  hybrid  matrix  exists, 
but  its  lxl  principal  submatrices  are  both 
singular  and  so  a  1-fold  gyration  is  never 
possible.  Hence  only  two  corners  and  the  con¬ 
necting  diagonal  of  the  2-diraensional  hypercube 
exist.  The  degenerate  cube  is  depicted  in 
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Figure  37  .  The  impedance  and  admittance  matrices 
are  undefined.  Both  existing  corners  are  associ¬ 
ated  with  an  IPR  matrix  of  degree  0  with  a  null 
Itermitian  part. 

Example  4. 

Consider  the  null  matrix  / 0  0\,  IPR  of  degree  0. 


Oi 


If  this  is  regarded  as  an  impedance  matrix  then  it 
is  realizable  by  two  short  circuits.  No  gyrations 
are  possible  end  so  the  two-dimensional  cube  de¬ 
generates  into  a  single  point  at  (0,0). 

As  an  admittance  matrix  it  is  realizable  by 
two  open  circuits,  and  has  a  degenerate  cube  con¬ 
sisting  of  a  single  point  at  (1, 1) . 

As  a  hybrid  matrix  T^(Z)  it  is  realizable  as 
a  short  across  port  2  and  an  open  circuit  in  port  1. 
This  follows  from  the  fact  that  in  this  case 


As  the  hybrid  ^(Z)  it  is  likewise  realized  by  a 
short  and  open  circuit.  In  the  last  two  cases  as 
in  the  first  two  the  cube  degenerates  into  the 


single  points  (1,0)  and  (0,1)  respectively. 


CHAPTER  IV 


General  PR  inunittance  Matrix  Synthesis 

4.1  INTRODUCTION  In  this  final  chapter  the 
gyration  operator  is  used  as  the  basis  for  a 
general  synthesis  procedure  for  PR  immittance 
matrices.  The  synthesis  procedure  derived  in 
Chapter  ii  is  a  special  case  of  this  general  pro¬ 
cedure.  In  this  chapter  the  symmetry  condition 
is  not  however,  required  and  both  non- symmetric 
and  symmetric  matrices  can  be  realized  by  this 
general  procedure. 

4.2  COMMENTS  ON  THE  DEFINITION  OF  A  PR  MATRIX 

Let  F(s)  be  a  matrix  function  of  thecomplex 

variable  s  all  of  whose  elements  f . .  are  rational 

ID 

with  real  coefficients. 

F  can  be  split  into  the  sum  of  two  matrices 

F  —  F  +  F 

H  SH 

where  F„  is  hermitian  and  F_„  is  skew-hermitian. 
n  SH 

This  follows  immediately  from  the  identity 
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F  +  F  *  F  -  F  * 
F  =  - : - —  +  - — 


Letting 


F  = 
H 


F  +  Ft* 


we  see  that  F  =  F  *  and  is  thus  hermit ian. 


Similarly  if 


F  -  F  * 
_ T 


then  F_„  =  -  F_  *  and  is  thus  skew-hermit ian. 

SH 

Condition  II  of  the  definition  of  the  PR  property 
for  matrices  requires  that  Re  x*Fx|  ^  0  for 
Re  h  -  0  and  arbitrary  ST.  Now 


x*Fxi  =  x*F„xJ  +  x*F  xl 
n  ori 


a  +  jp 


where  a  and  jp  are  real  and  pure  imaginary  since 
they  are  hermitian  and  skew-hermit  Lan  forms  re¬ 
spectively.  Hence  II  may  be  restated  as 
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II 


1 


The  hermitian  part  of  F  is  non-negative 
definite  for  Re  s  -  0, 


Ila,  lib  and  lie  can  be  restated  as 


F„  is  non-negative  definite  for  Re  s  s  0 
H 

f  .  has  no  poles  for  Re  s  >  0 

X  J 

i,  j  «  1,  2,  .  .  .  n. 

For  Re  s  =  0/f^j  has  simple  poles  and 
the  matrix  of  residues  is  hermitian  non¬ 
negative  definite. 


Likewise  IV  is  equivalent  to: 

IV*  A  positive  real  matrix  F  is  said  to 
be  IPR  if  on  the  line  Re  s  =  0,  F  is 
everywhere  skew-hermitian. 

Thus  the  rank  of  the  hermitian  part  of  an  IPR 
matrix  is  zero  everywhere  on  the  j-axis. 

As  will  be  seen  later,  this  restatement  of 
the  PR  property  in  matrix  terms  rather  than  the 
scalar  terms  previously  used  will  be  advantageous 
under  certain  conditions. 
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4.3  SYNTHESIS  PROCEDURE 

In  what  follows  we  develop  a  syntheses  pro¬ 
cedure,  using  the  r  operator,  </hich  is  capable 
of  handling  any  PR  immittance  matrix.  As  will  be 
shown,  any  such  matrix  falls  into  precisely  one 
of  two  classes  and  accordingly  the  procedure  will 
have  two  variants,  termed  Case  A  and  Case  B.  case 
A  will  be  shown  to  be  a  generalization  of  the  pro¬ 
cedure  already  given  in  Chapter  II  for  symmetric 
matrices  and  will  handle  all  symmetric  PR  as 
well  as  some  non-symmetric  PR  matrices,  case  B 
will  cover  all  remaining  PR  matrices. 

Suppose  we  are  given  a  PR  impedance  matrix 
£  which  is  to  be  realized.  As  a  first  step,  all 
imaginary  axis  poles  are  removed  by  splitting  the 
matrix  into  the  sum  of  pole  matrices  and  a  re¬ 
mainder  which  has  no  j-axis  poles.  All  of  these 
matrices  are  PR  and  the  pole  matrices  being  IPR, 
can  be  readily  synthesized  by  inspection  using 
inductors,  capacitors  and  possibly  gyrators  to¬ 
gether  with  ideal  transformers. 

A 

Assume  now  a  PR  impedance  matrix  Z  without 
j-axis  poles.  As  discussed  in  section  2.3,  by 
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removal  of  resistance  from  each  port  or  from  just 
one  port  a  situation  is  reached  where  the  herm¬ 
it  ian  part  Z  becomes  singular  at  some  point  on 
H 

the  j-axis  say  jo>0.  In  Lemma  5  below  it  is  shown 

that  there  exists  a  congruence  transformation 

which,  at  s  =  jco0,  places  a  zero  in  the  1, 1 

position  of  Z  or  else  places  a  singular  2x2 
H 

submatrix  in  the  first  principal  position  of  z  . 

H 

Lemma  5.  Let  F  be  an  n  x  n  matrix  which  is  not 
identically  singular  and  has  no  j-axis  poles. 

Let  F  be  PR  and  let  FH(j^o)  be  singular. 

Then  there  exists  a  real,  constant  non¬ 
singular  matrix  D  such  that  either 

i)  DmF„(jo)_)D  has  a  zero  in  the  1,1 
T  H  O  1  1 

posicion  (case  A) . 
or 

ii)  D^F^  ( jcoo)  D  has  a  2x2  submatrix  in  the 
first  principal  position  which  is 
singular  (Case  B) . 

Proof.  fh^coo^  s^n9ular  by  assumption.  Two 
cases  arise.  Either  Re  FH(lkoo)  is  singular  (Case  A) 
or  it  is  not  (case  B) . 


I 
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Case  A  and  Re  are  singular  at  jo)Q.  There 

exists  D,  real  and  non-singular  such  that 

DTRe  PH^<do^D  is  dia<?°nal  and  has  a  zero  in  the 

1,1  position.  But  F„  =  Re  F„  +  j  Im  FtI  and  Im  FTI 

«  H  H  H 

is  skew- symmetric,  a  property  which  is  preserved 
under  congruence. 

Hence  DTFH(ju)Q)D  has  a  zero  in  the  1,1  posi- 
t ion .  This  proves  case  A. 

Case  B  Re  FH(jo)Q)  is  non-singular.  By  the  PR 
property  of  F,  F„(jas  )  is  non-negative  definite 
(Condition  II1,  section  4.2).  Hence 

x*FH(jo>o)x|  =  3c* Re  FH(jo)o)xl  +  jx*im  F  ( joiQ)xl  *  o 
for  all  x. 

In  particular,  for  all  real  x  since  Im  FH(ja>0) 
is  akew-symmetric, 

xlm  F„(jm  )xl  =  0 
n.  o 

giving 

xRe  FH(ja)Q)xf  2s  o 

Hence  Re  FH ( j^0 )  is  non-negative  definite  and 
since  Re  FH(jcoo)  is  non-singular  by  assumption,  it 
is  in  fact  positive  definite.  Thus  there  exists 
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a  matrix  Q  which  is  real  and  non-singular  such 
that 


Let 


QTRe  s  1 

QTIm  FH(jcDQ)Q  =  N,  say 


where  N  is  real  skew-symmetric.  It  is  shown  in 
Appendix  2  that  there  exists  an  orthogonal  matrix 
V  such  that  if  the  real  skew-symmetric  matrix  N 
is  of  even  order  then  VTN  V  = 


and  if  N  is  of  odd  order,  then  V^N  V  = 
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In  either  case  some  of  the  may  be  zero.  Taking 
D  =  QV,  ve  obtain  = 


T 


r 


7S* — 'JST 


nr 


(assuming  that  n  is  even) ,  which  has  ones  in 
every  diagonal  position. 

But  F„(jaj_)  is  singular.  Hence  for  some  k, 

H  O 


is  singular.  Without  loss  of  generality,  let  this 
be  true  for  k  =  1.  Then  DmF„(jm  )D  = 

'  1  xl  O 


which  proves  case  B,  and  completes  the  proof  of 
Lemma  5. 
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f 

! 


Thus  after  Z„  has  been  made  singular  at  jcs 
H  o 

by  the  removal  of  resistance  the  matrix  is  found 
to  fall  into  Class  A  (all  those  for  which  Re  Z 

H 

is  also  singular  at  joo0)  or  Class  B  (those  for 
which  Re  is  positive  definite) .  We  observe 

Xl 

that  all  symmetric  matrices  immediately  fall  into 
Class  A  since  for  them  the  hermitian  part  is  real. 

We  now  apply  the  appropriate  congruence 
transformation,  as  given  by  Lemma  5  and  we  can 
thus  assume  without  loss  of  generality  that 

A 

DTZ  D  =  Z 

has  either  a  .hermitian  patt  whose  1, 1  element 
is  zero  or  else  it  has  a  hermitian  part  whose 
first  principal  2x2  submatrix  is  singular  at 

jcv 

This  congruence  transformation  by  D,  can  be 
compensated  for  by  also  applying  the  inverse 
congruence  transformation,  the  latter  being  real¬ 
ized  by  the  appropriate  use  of  ideal  transformers 
as  discussed  in  [6]. 

Thus  far,  the  synthesis  procedure  has  re^ 
sulted  in  the  network  of  Figure  12 


T 
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4.4  CASE  A  SYNTHESIS 

Assume  Z  to  be  an  n  x  n  PR  impedance  matrix, 

without  j-axis  poles,  with  z  (jm  )  = 

H11  ° 

If  z^  is  identically  zero,  we  have  the  following 
Lemma  6.  Let  F  be  an  n  x  n  PR  matrix  with  f^  =  0 
identically.  Then 

/f . ,  !  Ft  A  /  0  !  a  \ 


F21* 


0 

i  a 

i 

1 

1 

ai 

'  “1  "" "  - - 

1 

i  p22 

where  a  is  a  real  constant  row-vector,  and  F^2 
is  PR. 

Proof.  F22  is  PR  since  it  is  a  principal  submatrix 
of  a  PR  matrix.  F^  is  hermitian  non-negative 
definite  for  Re  s  -  0.  Since  f.^  =  0,  we  must 
have 


F12  +  (F21^t  =  ^  for  Re  s  ~ 

Let  s  be  real  and  positive.  Then  (F2l^* 
Hence  F12  +  (F2iOT  =  0  for  s  >  0,  and  so 

*12  =  ~(F21*  }T  f°r  a11  S' 

We  then  have,  F^2  -  F^|  =  for  Re  s  *  0. 


F21* 
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Let  s  be  on  the  j-axis. 

Then  f12  -  rxJ  =  "0  implies  that  F^  ra3l  on 
the  entire  j-axis. 

Hence  F^  =  a  a  real  constant  n-1  vector. 

This  completes  the  proof  of  the  lemma. 

If  z^  is  identically  zero,  we  may  by  Lemma 
6,  split  Z  into 


The  first  of  these  two  matrices  can  be  realized 


by  a  gyrator  and  ideal  transformers  (see  later) 
and  we  then  resume  synthesis  on  Z22*  Thus  we  may 
assume  without  loss  of  generality  that  z^  is 
not  identically  zero. 

Suppose  iiext  that  is  identically  zero 

H11 

on  the  j-axis.  Then  z^  is  IPR,  and  since  it 

cannot  have  poles  on  thej-axis,  (these  have  all 

been  removed),  z-j^  must  be  identically  zero.  But 

this  has  already  been  ruled  out.  Thus  we  may 

assume  that  z„  is  not  identically  zero  on  the 
H11 

j-axis. 
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Suppose  finally  that  Z  or  one  of  its  princi¬ 
pal  r inors  is  identically  singular.  We  then  have 
Lemma  7 .  Let  E(s)  be  £  PR  matrix  which  is  identi- 
cally  singular.  Then  F  is  congruent  to 


where  w  is  a  PR  matrix  and  a  is  a  real,  constant 
row-vector. 

Proof.  Take  s  -  1.  Let  F(l)  =  M  +  N  where  M  is 
real  symmetric  and  N  real  skew-symmetric.  Let 
C  be  a  real  constant,  non-singular  matrix. 


Then  CTF(1)C  =  CTM  C. 

But  C^F ( 1 ) C  is  singular. 

C  so  that  . 

/  0 


ctm  c 


Hence  so  is  M. 


Choose 


Then  [C^CsJc^^  has  a  zero  in  Re  s  >  0.  But 
this  implies  that  [ CTF  (s) C] ^  =  0  identically 
since  it  is  PR. 

Hence  by  Lemma  6 
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CtF(s)C 


I 


a 


w 


where  w  is  PR  and  a  is  a  real  and  constant  row 
vector.  This  completes  the  proof  of  the  lemma. 
Corollary.  if  F  JLs  PR  and  identically  singular, 
then  Fh  is  identically  singular. 

Proof.  There  exists  a  non-singular  c  such  that 


CTF  C 


a 


w 


\ 

/ 


Hence 


!' 0  |  <t>  \ 

!  WH  / 


which  proves  the  assertion. 

Thus  we  may  assume  that  none  of  the  principal 
submatrices  of  Z  is  identically  singular,  since 
if  this  were  the  case,  we  can  apply  a  congruence 
transformation,  remove  a  gyrator  section  and 
resume  synthesis  of  the  remainder. 


f 
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Then  ( jtoQ)  =  0  but  is  not  identically  zero. 

Z  has  no. j-axis  poles/  it  follows  that  Z  has 
only  the  j-axis  pole  possibly  due  to  b.  If 
0  a)  -c  oo  ,  then  Z  has  at  most  one  j-axis  pole 
in  only  z ^  and  this  pole  is  either  at  0  or  oo , 
depending  on  b.  On  the  other  hand  if  cdq  =  0  or 
oo  then  b  =  0  and  so  Z  then  has  no  j-axis  poles. 
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We  are  now  in  a  position  to  reduce  the  degree 
of  Z  by  the  removal  of  a  lossless  section. 


4.5  STATEMENT  AND  PROOF  OF  CASE  A  THEOREMS 
Theorem  18.  Let  Z(s)  be  an  n  x  n  PR 
impedance  matrix, 

a)  with  no  identically  singular  principal 
submatrices, 

b)  without  j-axis  poles, 

c)  with  2..  not  identically  zero  on  the 

H11  ~ 

j-axis,  and 

d)  z  =  0  at  jco  . 

H11  ° 


Let  Z 


!  j\ 


:  <t>/ 


e  the  IPR  scalar  b  is 


so  chosen  that  =  Then  Z  may  be  de¬ 

composed  in  the  series  parallel  manner  of  Figure 


1  where 

1. 

Z' 

is  PR  and  6Z'  =  6Z  -  6Z" 

0 

Z" 

is  IPR  and  6Z"  is  1  or  2 

3. 

Z" 

contains  at  most  2  reactors  and  1 

qyrator  plus  ideal  transformers. 
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Since  z ^  has  zeros  at  +jo)  ,  P-^Z)  has  poles  there. 
T1(Z)  is  PR.  Hence  by  Theorem  15,  F^(Z)  can  be 
split  into  an  IPR  matrix  P,  formed  by  the  re¬ 
moval  of  the  poles  at  +jcoQ  from  r^(Z),  plus  a 
matrix  Q  =  r^(Z)  -  P  which  is  PR  without  poles 
at  ±jco0- 

As  will  be  shown  under  the  proof  of  Pro¬ 
position  3,  p  ^  is  not  identically  zero  and  so 
Z"  =  p^ (P)  exists. 

Suppose  is  identically  zero.  Since  p^ 

is  IPR,  it  has  a  real  part  which  is  identically 
zero  on  the  j-axis.  It  follows  from  the  corollary 
to  Theorem  9  that  z^  =  r^(Pn  +  also  has  a 

real  part  which  is  identically  zero  on  the  j-axis. 

But  this  contrary  to  assumption  c  in  the 
statement  of  this  theorem.  Hence  cannot  be 
identically  zero  and  so  Z'  =  T^(Q)  exists. 

By  Theorem  2,  both  Z'  and  Z"  are  PR.  Since 
P  and  0  share  no  poles  it  follows  from  Equation  1, 
that  6 (P  +  Q)  =  6P  +  6Q.  By  Theorem  1,  the  degree 
is  invariant  under  a  gyration,  so 


6Z  =  6 (P  +  Q)  =  6Z*  +  6Z". 


Hence 


6Z'  =  az  -  az" 
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We  shall  refer  to  these  as  propositions  1,  2, 
and  3  respectively. 

This  theorem  is  the  basis  of  the  Case  A  pro¬ 
cedure.  In  order  to  compensate  for  the  addition 
of  b(s)  we  must  also  add  -b(s)  to  z^i*  It  turns 
out,  analagously  to  the  scalar  case,  that  this 
subsequent  negative  element  can  be  incorporated 
into  a  perfectly  coupled  transformer.  This  is 
proved  in 

Theorem  19.  Let  Z,  Z  and  b  be  as  defined  in 
Theorem  18.  Then  after  the  application  of  that 
theorem  to  split  Z ,  it  is  always  possible  to  in¬ 
corporate  the  compensating  negative  reactor  into 
a  perfectly  coupled  transformer. 

Proof  of  Theorem  18. 

Proposition  1. 


Form  r^Z)  =  (zi;L)~  /  1  \  (1 


-Z12> 
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This  completes  the  proof  of  Proposition  1. 
Prior  to  proceeding  with  the  proofs  of  the  re¬ 
maining  propositions  we  make  the  following  ob¬ 
servation. 


Let  G  be  a  real  skew-symmetric  n  x  n  matrix. 
Then  both  G  and  -  G  have  a  non-negative  definite 
hermitian  part  and  so  they  are  PR.  Thus  P  +  G 
and  Q  -  G  are  also  PR  and  it  is  thus  permissible 
to  split  F1(Z)  into  P  +  G  and  Q  -  G.  Since 
6G  =  0  it  follows  then  that  all  the  arguments 
given  in  the  above  proof  of  Proposition  1  are 
equally  valid  if  Z"  =  r-^P  +  G)  and  Z'  =  r^(Q  -  G). 

G  will  be  used  to  a  definite  advantage  in  the 
proof  of  Theorem  19  to  follow. 


First  let  0  <  a)  <  oo 

o 

Then  the  pole  matrix  removable  from  T^(Z)  is 
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(s  -  jo)Q)U 


(s  +  jcno)ul 


s=DaV 


s=-3a>. 


s  -  ICDr 


s  +  3m 


V 

s  -  jui 


V* 

S  +  jcD 


Now  u  is  of  rank  1.  Hence  the  residue  matrices 
V  and  V*  are  of  rank  1.  Thus  by  Equation  1,  P  is 
of  degree  2,  which  means  that  by  Theorem  1, 

Z"  =  is  of  degree  2.  P  is  IPR,  and  so  the 

rank  of  the  hermitian  part  of  P  is  zero  on  the 
j-axis  (Condition  iv\  section  4.2) 

Thus  by  Theorem  2,  Z"  is  PR  and  by  the  corol¬ 
lary  to  Theorem  9,  Z”  is  IPR.  Hence  Z"  is  of  de¬ 
gree  2  and  IPR. 

Let  o)Q  =  0  or  00 


Suppose  for  definiteness  that  a>  =0. 

o 

(The  case  o)Q  =  00  can  be  analyzed  as  the  dual) 


Then  P  = 


.  U  is  of  rank  1,  and  so 


P  is  of  degree  1.  P  is  IPR.  Hence  Z"  =  r^(P) 
is  of  degree  1  and  IPR. 
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Thus  for  0  —  oo  —  od  ,  Z"  is  IPR  and  of  degree  1  or 


This  completes  the  proof  of  Proposition  2. 
Before  proving  proposition  3,  we  prove 
Lemma  8.  Let  F  be  partitioned  as  follows 


fll  ! 

i 

i 

- 1. 


where  f  is  a  scalar. 


Let  F  be  PR  with  f„.  (ja>  )  =0,  but  not  identi- 
__  «n  o 


cally  zero.  Then 


=  or  +  j-pr 


F21(jwo)l  =  _aJ  +  3‘Pl 


where  o'  and  (3  are  real  (n  -  1)  vectors. 

Proof.  F..  is  non -negative  definite  on  s 
*  H 

since  f  (jo)  )  =  0,  we  have  at  ju)  , 

H11  °  ° 


ju>  and 


f tj  ® f tt  ® /  k  —  2,3,  .  .  .  n 

Hlk  Hkl 
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V  V* 

s  -  jo>0  s  +  ja>0 


Let  V  =  M  +  jN 


Then  P 


2s 


2  2 
s  +  CO 

o 


kl 

1 

1 

1 

M12 

— 

- 

1 

— 

iM2ll 

1 

1 

1 

M22 

2co, 


s  +  co. 


11 
- 1 


N21^ 


Let  V  = 


v  i 

11  i 

i 
i 
i 

- r 

i 
i 


v^2\ 


\ 


V2l‘ 


V 


22/ 


Then  vn  =  (B  -  jo,0)  - 


11 


=  7  say, 


s=jm 


which  is  real  and  positive  since  is  PR. 

i.e.  =  7  and  n^  =  0. 


V12  =  -(B  -  j„,o)  z 


12 

11 


s=3co. 


=  -7 {a  +  jf3)  by  Lenuna  8,  and  so 


Mi 2  =  -Va,  N12  =  -7(T. 


Similarly  v2jl  =  7  (— cx|  +  j^{  )  by  Lemma  8 
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Thus  Z"  has  been  dplit  into  the  three  matrices 
given  in  the  preceeding  equation.  The  first  of 
the  three  matrices  of  Z"  may  be  realized  as  shown 
in  Figure  13.  The  second  term  o\  Z"  can  be  real¬ 
ized  as  shown  in  Figure  14. 

The  third  part  of  Z"  namely 


(  0 
\“a{ 


\ 


/ 


can  be  realized  by  one  gyrator  and  a  congruence 
transformer.  This  follows  from  the  fact  that  is 
skew  symmetric,  real,  of  rank  2.  Hence  there 
exists  a  non-singular  matrix  V  (see  Appendix  2) 
such  that 


which  is  realized  in  Figure  15.  The  overall  net¬ 
work  for  Z"  is  shown  in  Figure  16. 


X  -£T 

4 


Xw 
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An  alternate  realization  for  Z"  is  possible. 
Let  9  be  a  real  n-1  vector.  Observe  that  both 


and  -G  are  real  skew-symmetric  and  hence  they  are 
PR.  Instead  of  splitting  T^(Z)  into  P  and  Q  we 
split  it  into  P  +  G  and  Q  -  G  as  was  anticipated 
in  the  comments  after  the  proof  of  Proposition  1. 


Let 
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Z"  of  the  alternate  realization,  has  been  split 
into  three  matrices.  The  first  of  these  terms 
is  a  capacitor  in  series  with  port  1  of  the  net¬ 
work  of  the  remaining  term3. 

The  second  term  is  a  single  perfectly  coupled 
transformer  shown  in  Figure  17. 

The  third  term  can  be  realized  by  a  single 
gyrator  and  a  congruence  transformer,  since  it 
is  congruent  to 


The  overall  network  for  Z"  is  given  in  Figure  18. 

If  oo  =  0  or  co 
o 


For  definiteness  assume  a>  =  0 


analyzed  as  the  dual) . 

Then  U  =  (z^)  ^  /l  \  (1 


(u>  =  oo  can  be 

o 


) 
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The  pole  matrix 


s 


V 

s 


say. 


v, 

X 


1 


=  Y  which  is  real 


and  positive. 


~SZ12 

Z11 


s=0 


But  Z12(0)  is  real.  Hence  by  Lemma  8,  Z12(0) 
and  Z21(0)|  =  -q|.  Thus 


a 


=  -7a 


V_ ,1  -  -Yal 

Z  a. 

Finally 

v22  =  Tala 

f1 

Hence  P  =  -  —  - 
s 

^-ai 
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Then  Z"  =  f^P) 


Thus  Z"  can  be  realized  by  a  single  inductor,  a 
single  gyrator  and  a  congruence  transformer,  as 
shown  in  Figure  19. 

Hence  Z"  contains  at  most  two  reactors,  one 
gyrator  and  ideal  transformers. 

This  proves  Proposition  3  and  completes 
the  proof  of  Theorem  13. 

Proof  of  Theorem  19. 

We  recall  that  Z  was  partitioned  as  follows 


where  z^  is  a  scalar. 

Then  since  z  (joj  )  =  0,  we  have 
H11  ° 


I 
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2ll(ja3o)  =  jT) 


where  tj  is  a  real  constant. 


Let  0  <  to  <  co 
o 


If  r)  >  0  let  B  = 


i  : 


i  4» 


If  7)  <  0  let  B  =  - 


<t>i-  !  <t> 


If  7)  =  0  let  B  =  <t>. 


If  cd  =  0  or  oo 
o 


Then  t)  =  0.  Let  B  =  <t>. 

in  each  of  the  above  cases  the  matrix  B  is  PR. 
Recall  that  Z  =  Z  +  B.  Then  Z  is  PR  and 


*11  (joy  =  0. 

Case  1  z  (jco  )  =  j r) 
H11  ° 


7)  >  0. 


Then  z„,  =  zni  +  -  ha3  a  zero  at  co  and  a 

11  11  s  o 

pole  at  the  origin. 


Recall  that  p  .  =  — 5 - 2  •  which  is  subtracted 

1  s  +  CD 

o 


from 


1 

Z11 


leaving 
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2 

"  ^  *0  .  S 

Z11  2s7  27 

and  utilizing  the  alternate  relization  for  Z" 
(Figure  18)'  we  obtain  the  network  of  Figure  20 
where  the  series  capacitor  -oi  t)  /  s  i-s  the  com¬ 
pensation  for  adding 


to  Z.  Now  the  three  capacitors  of  Figure  20  can 
be  combined  into  a  three-terminal  network 
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which  has  a  pole  at  co  given  by 

-an 

%  +  27rj 

Since  this  term  is  Pr  and  ?i  <0,  we  have  that 
+  2 7rj  >  0.  From  Z"  we  obtai:  the  pole  at  oo 
of  value  s/,27  (see  the  proof  of  Proposition  3 
of  Theorem  18) .  The  three  terminal  matrix  T  is 
now: 


which  is  PR  since  ao  +  27-n  >  0  and  0  <cco  <oo  . 

o  '  o 
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Again  T  is  of  degree  1  since  it  has  rank  1.  T  is 
realized  by  the  three  terminal  network  of  Figure 
23  which  is  a  single  inductive  element,  and  Z  is 
realized  by  Figure  24.  This  completes  the  proof 
of  Theorem  19. 

Prior  to  commencing  with  a  discussion  of  Case 
B,  we  make  the  following  observations. 

1)  Referring  to  Figures  22  and  24,  either 


or 


Z' 


/i 


<i> 


\ 


+ 


s  n 

aS  +  277) 
o 


now  remain  to  be  synthesized.  If  coQ  =  0  or  co , 
then  Z'  remains.  The  cycle  is  now  resumed  by 
splitting  off  any  j-axis  poles  followed  by  a  re- 

,  A 

moval  of  resistance  until  a  new  Z  is  produced 

A 

for  which  Z  is  singular  at  some  point  on  the 

ri 

j-axis.  Either  a  Case  A  or  a  Case  B  cycle  now 


commences . 
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2)  We  note  that  a  buffer  of  degree  1  or  2 
was  removed  4 rom  Z  resulting  in  a  degree  reduction 
of  1  or  2.  Thus  the  Case  A  cycle  is  minimal  in 
the  sense  of  Tellegen' s  Minimal  Theorem  stated  in 
section  2.3.  As  will  be  shown  later,  the  Case  B 
cycle  reduces  the  degree  of  Z  by  2  with  the  re¬ 
moval  of  a  buffer  of  degree  2  and  so  it  is  also 
minimal.  Thus  a  given  matrix  of  degree  d  will  be 
synthesized,  by  this  method,  by  a  network  contain¬ 
ing  precisely  d'  reactive  elements. 

3)  We  note  that  at  most  1  gyrator  appears 
in  the  buffer  for  a  Case  A  cycle,  whether  the  de¬ 
gree  of  the  buffer  is  1  or  2.  We  wil>l  show  that 
at  most  one  gyrator  is  required  per  unit  degree 
reduction. 

4.6  CASE  B  SYNTHESIS 

Assume  that 


/* ii  ! 

1 

i 

?12\ 

1 

\~Z21  i 

-22/ 
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is  an  n  x  n  PR  impedance  matrix.  is  2  x  2. 

Let  Z  have  no  j-axis  poles,  and  assume  as  in  Case 

/Si 

A  that  none  of  the  principal  submatrices  of  Z  is 
identically  singular.  Since  Z  falls  under  Case 
B,  then  Z„  (jco  )  is  singular,  but  Re  ZTT  (jco  ) 

Hn  °  Hn  ° 

is  positive  definite.  This  immediately  means 

ssi 

that  cd  /  0  or  oo  since  at  those  points  Z  is  real 
o 

and  so  Z  (jco  )  is  real  and  singular,  which  im- 
H11  ° 

plies  that  a  Case  A  synthesis  is  possible.  Sim¬ 
ilarly  Z  cannot  be  identically  singular  on  the 
hll 

j-axis  since  then  it  is  singular  at  coo  =  0  and  oo . 

A  case  A  synthesis  is  then  also  possible. 

Thus  we  may  assume  without  loss  of  generality  that 

a'  #  m  # 

for  Case  B,  Z  is  an  n  x  n  PR  impedance  matrix, 
without  j-axis  poles,  with  no  identically  singular 

fSj 

principal  submatrices,  Z„  not  identically  sing- 

11 

/V 

ular  on  the  j-axis,  Z  (jco  )  singular,  Re  Z  (jco  ) 

W11  °  “ll  ° 

not  singular  and  0  <  cdq  <  oo  . 

We  now  follow  a  procedure  similar  to  the 
Brune  method,  and  add  to  Z^  a  2  x  2  matrix  B(s) 
such  that  Z^  +  B  is  singular  at  jcoQ .  The  exist¬ 
ence  of  such  a  matrix  is  guaranteed  by 
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Lemma  9 .  Lot  Z^  be  a  2  x  2  PR  matrix  with 
Zn(j®o>  =  C  +  jD  -  j  (E  +  jF)  where  C  and  E 
are  real  symmetric,  D  and  F  are  real  skew- 
symmetric  and  where  C  +  jD  is  positive  semi- 
definite  (i.e.  singular)  but  C  is  non-singular. 

Then  there  exists  a  2  x  2  PR  matrix  B ( s )  such 

that 

i)  zn^a)0)  +  B  ( joDQ)  is  singular 
ii)  B(s)  can  be  realized  by  at  most  1  in¬ 
ductor,  i  gyrator  and  1  ideal  transformer. 

Proof.  Let  B(jtto)  =  jA(C  +  jD)  +  j  (E  +  jF)  where 
A  is  as  yet  unspecified.  Then 

Zn(joOo)  +  B(ju>o)  =  (1  +  jA)  (C  -f  jD) 

which  is  singular  as  required  by  i)  of  the  lemma. 
To  give  B(s)  the  properties  required  by  ii)  of  the 
lemma,  A  is  chosen  as  follows.  C  is  non-singular 

A/ 

by  assumption.  Since  is  PR,  C  is  also  non¬ 

negative  definite.  Hence  C  is  in  fact  positive 
definite.  Thus  there  exists  a  real,  non-singular 


N  such  that 


I 
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C  =  NtN  and 


E  =  KL 


0 


r 


Then 


AC  +  E 


Suppose  that  p  >  q 


Take  A  =  -q.  Then  AC  +  E 


0 


N 


0 


/ 


which  is  positive  semi-definite. 
Take  B(s)  =  (AC  +  E)  -  (AD  +  F) . 


o 


Then 


B(jo)Q)  =  jA(C  +  jD)  +  j(E  +  jF) 

as  required  by  i) .  Moreover  B(s)  is  realizable 
by  a  gyrator,  a  single  reactor,  (since  XC  +  E 
is  positive  semi-definite)  and  ideal  transformers 
as  shown  in  Figure  25.  Thus  ii)  is  satisfied. 
Suppose  p  <  q 

Take  A  =  -p  and  follow  the  above  argument.  We 
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again  obtain  a  B(s)  which  satisfies  i)  and  ii) . 
Suppose  p  =  q. 

Take  X  =  -  p.  Then  XC  +  E  is  a  null  matrix.  Let 
B  (s)  =  -  (XD  +  F)  . 


Hence  B(joo^)  -  jX(C  +  jD)  +  j  (E  +  j F )  as  re¬ 
quired  and  B(s)  is  realizable  by  a  single  gyrator. 

Thus  in  every  case  B(s)  exists  which  satisfies 
i)  and  ii)  of  the  lemma. 

This  completes  the  proof . 

We  conclude  from  this  lemma  that  B(s),  a 
2  x  2  PR  matrix,  can  be  found,  whose  addition 
makes  )  singular.  This  is  an  extension 

of  the  Brune  procedure,  and  following  the  Brune 
tradition,  B(s)  will  be  compensated  for  at  a  later 
stage  by  the  addition  of  -B(s),  which,  as  will  be 
shown,  can  be  incorporated  into  a  Pr  network. 

Let 


Z 


Z  + 


\ 


B 


♦  \ 
♦/ 


Then  Zn  { ju)Q)  is 
cally  singular. 


singular  but  Z^  is  not  identi- 

'V 

Since  Z  has  no  g-axis  poles, 


Z  has  only  the  j-axis  poles  contributed  by  B, 
and  these  occur  only  in  Z^.  Since 

B  =  ^-{AC  +  E)  -  (AD  +  F),  if  AC  +  E  is  not  the 
o 

null  matrix,  Z^  has  a  rank  1  pole  at  oo,  and  if 
AC  +  E  is  null,  then  Z^  has  no  poles.  We  are  now 
in  a  position  to  reduce  the  degree  of  Z  by  the 
removal  of  a  lossless  section. 

4.7  STATEMENT  AND  PROOFS  OF  CASE  B  THEOREMS 

Theorem  20.  Let  Z  be  an  n  x  n  PR  impedance 
matrix 

a)  without  j-axis  poles 

b )  without  identically  singular  principal 
submairices 

c)  Z„  not  identically  singular  on  the 

H11 

entire  j-axis 

A/ 

d)  Z„  (j(JD  )  singular 

H11  °  - - 

e)  Re  ZTT  (jeo  )  not  singular  and 

H11  ° 

f)  0  <  tt>Q  <  oo  . 
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where  B  is  IPR,  2x2  and 


so  chosen  by  Lemma  9,  that  is  singular. 

Then  Z  may  be  decomposed  in  the  manner  of  Figure 
26  where 

1)  Z*  is  PR  and_  6Z*  =  6Z  -  2 

2)  Z"  is  IPR  an<d  6Z"  =  2 

3 )  Z "  consists  of  an  inductive  portion  of 
degree  2,  at  most  one  gyrator  plus  ideal 
transformers . 

We  will  refer  to  these  as  Propositicns  1,  2,  and 
3  respectively. 

This  theorem  is  the  basis  of  the  Case  B  pro¬ 
cedure.  In  order  to  compensate  for  the  addition 
of  B(s)  we  must  also  add  -B(s)  to  Z.^.  It  turns 
out,  analagously  to  the  scalar  case,  that  this 
subsequent  negative  matrix  can  be  incorporated 
into  a  PR  6-terminal  network.  This  is  proved  in 

A/ 

Theorem  21.  Given  Z,  Z  and  B  as  defined  in 

Theorem  20  and  Lemma  9.  Then,  after  the  application 
of  Theorem  20  to  split  Z,  it  is  always  possible  to 
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incorporate  the  compensating  -B(s)  into  a  PR  net¬ 
work  .  Proof  of  Theorems  20  and  21  now  follow. 
Proof  of  Theorem  20. 

Proposition  1 


Form  ri  2 (Z) 


) 


lz 


21 


I* 

!  *  \ 

+ 

1 

r 

!  Z22/ 

Since  is  singular  at  +ja>Q,  T ^  ^{7.)  has 

poles  there.  By  Theorem  2,  r ^  2(Z)  PR' 

Hence  by  Theorem  15,  T,  9(Z)  can  be  decomposed  in- 
to  an  IPR  matrix  P,  formed  by  the  removal  of  the 
poles  at  ju>o  from  2^'  an'^  a  nlatr^x 
Q  =  f  ~(Z)  -  P  which  is  PR  without  poles  at 


+jU>^. 


Partition  P  into  /p..  1 


11  j  12 


\P21  !  P22/ 


where  P 


11 
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is  2  x  2.  Do  likewise  for  Q.  As  is  shown  below 
in  the  proof  of  Proposition  3,  P^  is  not  identi¬ 
cally  singular  and  so  Z"  =  ^ (P)  exists. 

cannot  be  identically  singular,  for  suppose  it  is. 

By  the  corollary  to  Lemma  7,  Q  is  then  identi- 

H11 

cally  singular.  Since  P^  is  IPR  it  follows  that 

% 

the  hermitian  part  of  (P.^  +  is  singular  on 

the  entire  j-axis.  Hence  by  the  corollary  to 

Theorem  9,  =  T  ^  2^11  +  ^11^  ^as  3  hermitian 

part  which  is  singular  on  the  entire  j-axis.  This 
contradicts  assumption  c)  of  this  theorem.  Hence 
cannot  be  identically  singular  and  so 
Z  =  F  ^  2 (Q)  exists . 

By  Theorem  2,  both  Z'  and  Z"  are  PR.  P  and 
Q  share  no  common  poles.  Hence  it  follows  by 
Equation  1  that 

6 (P  +  Q)  =  6P  +  6Q 

But  by  Theorem  1,  the  degree  is  invariant  under  a 
gyration,  so 

6Z  =  5P  +  5Q  =  5Z'  +  5Z" 


It  will  be  shown  below,  in  the  proof  of  Proposition 


-115- 


3,  that  6Z"  =  2.  Hence 


5Z '  =  5Z  -  2 


This  completes  the  proof  of  Proposition  1. 

We  make  the  following  observation.  Let  G  be 

a  real  skew-symmetric  matrix.  Then  both  G  and  -  G 

are  PR.  If  instead  of  splitting  T,  0(Z)  into  P 

1. 1  z 

and  Q,  we  split  it  into  P  +  G  and  Q  -  G,  then 
both  P  t  G  and  Q  -  G  are  PR. 

If  we  choose  G  so  that  it  has  the  form 


=  u 


4> 


( 

0  I  F 


where  Fis  (n-2)  x  (n-2)  then 


by  the  same  arguments  as  in  the  proof  of  Prop¬ 
osition  1,  P^  and  Q  are  not  identically  sing¬ 
ular.  Hence  Z"  =  T,  (P  +  G)  and  Z'  =  F ,  „ (Q  -  G) 

x  t  Z  J- 1  £ 

will  be  defined,  PR  and  of  the  same  degrees  as 
in  Proposition  1.  G  will  be  useful  in  Proposition  3. 
Prior  to  proving  Proposition  2,  we  prove 
Lemma  11  Let  M  be  a  2  x  2  matrix.  Let  A  and 


B  be  2  x  r  matrices.  Then 
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N  =  / \  M  ( I ^  J  B)  has  the  same  rank  as  M. 


Moreover  if  A  =  B  and  M  is  hermitian  positive 
semi-definite,  then  N  is  hermitian  positive  semi- 
definite  . 


2-element  row  vectors.  Then 


K  say  where  the  a\  are  linear 
1  l 


combinations  of  m^  and  m^ .  Let  K  =  (k^J  k2*  ) 
where  k^i  and  k,^  are  column  vectors  of  order  r+2. 
Then  K  (I  :  B)  =  (k^  k  |  .  ,  ,b  |  )  =  N 

where  the  b^  are  linear  combinations  of  k^  and 
.  Clearly  the  rank  *>f  N  equals  the  rank  of  K 
and  the  rank  of  K  equals  the  rank  of  M. 
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This  proves  the  first  part  of  the  lemma.  If  A  =  B, 


let  Q 


=  /• 


l 

I 

i 

f  • 
I 

i 


\  .  Then  Q  is  non-singular. 


Hence  N  is  congruent  to 


Qt*N  Q 


M 


I 

» 

i 

- 1  - 
i 
I 

i 


which  completes  the  proof  of  the  lemma. 

Proof  of  Proposition  2 

Recall  that  0  <  0)  <  oo 

o 


Define  U  =  /  I~  \ 


\Z2V 


-1 

'll 


(I. 


-Z12} 


Zu(jo>0)  a  singular  2x2  matrix.  Hence  its 

rank  is  zero  or  1 .  If  its  rank  is  zero,  then 

zll^coo)  =  ®  and  a  Case  A  synthesis  is  possible. 

Thus  we  may  assume  that  Z.^(ja)o)  ran^ 

at  s  =  jo)  is 
o 


The  residue  matrix 
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IT 


•#r7^^r 

% 


V  =  (s  -  j^0)U 


s=Jtor 


=  3  ~  jmo  I1 2  \adj  211  (I2  :  -Z12> 

det  Z,  , 


s=jo> 


Since  Z^  is  2  x  2  rank  1,  adj  Z  ^  has  rank  1.  It 
thus  follows  by  Lemma  11  that  V  is  of  rank  1. 

The  pole  matrix  removable  from  ^{Z)  is: 


(13) 


s  -  ju> 


s  +  jco 


which  is  thus  of  degree  2.  Moreover  P  is  IPR 
by  Theorem  15.  It  therefore  follows,  by  Theorems 
1  and  2  and  the  corollary  to  Theorem  9  that 
Z"  =  2(P)  is  IPR  and  of  degree  2. 

This  proves  Proposition  2. 

Proposition  3. 

The  residue  matrix  at  s  =  jWQ  in  f ^  2(Z)  is 


(s  -  j«>0) 
det  Z,, 


I2  \adj  Zu  (I2  : 


■Z12> 


s 
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which  is  rank  1,  as  was  noted  above  in  the  proof 

of  Proposition  2.  But  V  is  the  residue  matrix 

at  a  j-axis  pole  of  a  PR  matrix.  Hence  V  is 

hermitian  non-negative  definite  by  condition 

11^  ,  (section  4.2).  Partition  V  into 
c 


i 


I 


\ 


/ 


where  is  2  x  2.  Then  is  also  hermitian 

non-negative  definite.  Since 


V 


11 


(s  “  j<0 

det'Zll  adi  Z11 


S=ju) 


o 


and  since  is  2  x  2  rank  1,  it  follows  that  V 

is  in  fact  hermitian  positive  semi-definite  and 
must  therefore  be  of  the  form 


V 


11 


11 
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where  a^  >  0.  Moreover  b^  f-  0  since  if 

~  0/  then  z^^0u>o)  is  symmetric  and  singular 
and  so  a  Case  A  synthesis  is  possible. 

Having  established  the  form  which  must 

have  we  now  determine  the  form  of  V.  Let  C  and 
D  be  real  2  x  n-  2  matrices.  If  we  take 


then  V  is  hermitian,  and  by  Lemma  11,  it  is  posi¬ 
tive  semi-definite,  rank  1  as  required.  Let 
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Recall  that  P(s)  = 


Let  V  =  M  +  jN 


V 


V* 


S  -  jcD  S  +  jCD 

J  o  J  o 


(see  equation  13) 


Then  P  = 


2sM  -  2cd  N 
_  o 

2  2 

o 


Letting 


ai.d 


we  obtain 


/ 


sa 


K  = 


11 


sa,  „  -  a) 

12  o  12 


bi  i  \ 


a122  +  b122 

ySa12+“ob12  g<~  au  ) 


1^ 


L  = 


11 


CD  a  +  sb,  n 
o  12  l2 


\ 


U)oa12  ~  o  ( 


a  2  .  K  2 
12  +  b12 


12  o '  a. 


K 


KC  -  LD 


(14)  P  = 


CmK 


+dtl  I 


ctk  c  +  dtk  D 


“  ctl  d  +  V  c 


2  ,  2 
s  +0) 

o 


Hb 


We  now  add  to  P  the  matrix 


where  P  is  real,  skew-symmetric  (n  -  2)  x- (n  -  2) . 
As  noted  in  the  comments  after  the  proof  of  Prop¬ 
osition  1,  the  addition  of  this  matrix  to  P  can 
be  offset  by  subtracting  it  from  Q.  Then 


which  yields 


(15) 


25"  = 


where 


M 


ksostl'Mekw**,-*1 


-123- 


15  a  real  2  x  n-2  matrix 


-124- 


The  first  part  cf  Z"  consists  of  a  non-singular 
inductive  matrix  of  degree  2  since  a^  >  0  and 
£>12  /  0.  Hence  this  section  of  Z"  can  be  real¬ 

ized  by  2  inductors  plus  a  congruence  transformer. 

By  the  appropriate  choice  of  F>  the  second  part 
of  Z"  can  be  realized  by  1  gyrator  plus  a  congrui- 
ence  transformer.  To  see  this,  suppose  first 
that  M  =  is  of  rank  2  and  without  loss  of  gener- 

W 

03 

ality  let  —  =  1.  Let  F  =  m^H^  -  m2|m^. 

Observe  that  F  is  a-  skew  symmetric  real 
(n  -  2)  x  (n  -  2)  matrix  as  required.  Then  the  con¬ 
stant  part  of  Z"  is  congruent  to 


4> 
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where 


I 


0 


m. 


\ 


H 


0 


V 


1  I  -in, 
l  1 

-  -  + - 


/ 


Thus  the  constant  part  of  Z"  can  be  realized  by  1 
gyrator  and  ideal  transformers.  On  the  other  hand, 
if  M  is  of  rank  1,  i.e.  M  =  /  m"..  \  then  letting  F 

Vs] 

be  the  null  matrix  yields  the  result  that  only  1 
gyrator  plus  a  congruence  transformer  is  required. 
This  completes  Proposition  3  and  Theorem  20 . 

The  connections  of  Z'  and  Z"  will  be  discussed  in 
the  proof  of  Theorem  21  which  now  follows. 

Proof  of  Theorem  21 

As  required  for  a  Case  B  synthesis  Z  (jco  )  was 

H11  ° 

a  singular  2x2  matrix  where  0  <  ■<  oo  . 

We  then  added  B(s),  a  2  x  2  matrix,  chosen 
according  to  Lemma  9,  such  that 


B  (s)  = 


0) 


( \C  +  E)  -  (ID  +  F) 
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has  a  real  symmetric  part  XC  +  E  of  rank  1,  a 
real  skew-symmetric  part  XD  +  ^  and  has  the 
property  that  2’^  +  B  is  singular  at  s  -  jo)o. 


Now  r,  AZ  +1-  —  +-- 


=  (Z,,  +  B) 


Hence  Z^  =[(Z.^  +  B)  ^  -  P^  ]  ^  where  P  is  the 


pole  matrix  in  2 


Z  + - \ - 


at  jo>  (see 
o 


Equation  14) .  i.e. 


Z11  =  <Z11  +  BM  I  "  pn^zii  +  B)  3 


From  Bquation  14,  we  see  that  near  s  =  oo, 

P,  •,  4*  ~  ,  where 
11  s 


2  ,  v  2 

a12  +  b12 


is  positive  definite  since  >  0  and  b12  ^  °* 

Moreover,  Z  (and  hence  Z^)  contains  no  j-axis 
poles  by  assumption.  Tlsus  near  s  =  oo , 

Pll(Zn  +  B)  =?  A(XC  +  E) 

"L  JL"  o 

and  it  follows  then  that  Z^  has  a  2  x  2  pole  at 
given  by: 

M  =  ~(XC  +  E)  (I  -  A(XC  +  E)  f 1 
o  o 

which  can  be  removed  from  Z*  leaving  behind  a 
PR  remainder. 

Next,  by  Equation  15,  Z"  has  a  2  x  2  pole 
matrix  in  its  1,1  2,2  position 


which  can  be  removed  from  Z"  leaving  a  PR 


remainder, 
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Thus  the  connections  for  the  realization  of  Z 
are  as  shown  in  Figure  27 . 

The  addition  of  -B  is  to  compensate  for  the 
addition  of  B  to  Z.  Observe  the  common  terminals 
Ig  and  2g.  There  aie  in  fact  three  external 
terminals  (la,lb,lc)  associated  with  port  1. 
Similarly  for  port  2.  This  gives  a  total  of 
si:*  termina  Is  and  two  commons  for  the  network  made 
up  of  -B,  N  and  M. 

-B  is  not  PR,  but  it  can  be  absorbed  into  a 
6-terminal  (4  port)  matrix  given  by 
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shown  in  Figure  28. 

To  see  that  T  is  PR  we  proceed  as  follows, 


-B  +  N  =  [~  (XC  +  E)  4-  XD  +  F]  +  |  A"1 
o 


=  f  A-”1  [ I  -  |  A (XC  +  E)  ]  +  XD  +  F 


Define  Q  =  I  -  ^  A(XC  +  E) 


CD 


and  G  =  XD  +  F 

Then  -B  +  N  »  -|  A-1Q  4-  G 


Also,  N  +  M  =  |  A'1  +  |  (XC  +  E) [I 

o 


^  A (XC  +  E) ] 
o 


=  f  A*"1  [  I  -  |  A  (XC  +  E)  ]  1 
O'* 


=  §  A'V1 


Hence  T  = 


/  -1~ 

1 

-i  \ 

A  XQ 

1 

1 

i 

A  \ 

-1 

T 

i 

-1-1 

A 

i 

i 

A  Q  ^ 

\* 


4> 


s  -1  -3 

We  note  first  that  ^  A  Q  *  is  PR  since  it  was 


formed  from  the  addition  of  M  and  N,  both  PR. 
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* 


A  1  being  the  residue  at  a  j-axis  pole, 

must  be  hermitian  non-negative  definite  by  condi- 

1  s  -1  -1 

tion  IIIC.  By  condition  I,  A  Q  must  be 

real  for  s  real.  Hence  A  ^"Q  ^  is  symmetric  non¬ 
negative  definite. 

Since  both  Q  and  Q  ^  exist,  it  follows  that 
Q  ^  is  non-singular.  A  is  positive  definite. 

Hence  A  and  its  inverse  QA  are  both 

symmetric  positive-definite  2x2  matrices.  We 
then  have: 


_s 

2 


/  "1/s 

|  -1  \ 

A  Q 

1  A  \ 

-1 

t  -1  -1 

\A 

> 

a 

/ a-1/2  QA  a"1/2 


=  ■§  R 
2  T 


where 


\ 


\ 


!  A1/2  fll/2 


R 


/  .-1/2 


R  = 


♦  \ 


V 


!  A'1/2 
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But 


A"1/2  QA  a-1/2 

\  1 


\ 


aV2  a-V1  A1/2 


* 


\ 


A1//2  A_1Q_1  A1//2 


where 


Hence  the  reactive  part  of  T  is  congruent  to 


s 

2 


I 

i 

I 
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I 

I 

I 


4> 


A-V1 


\ 

/ 


and  can  therefore  be  synthesized  by  two  inductors 
and  a  congruence  transformer.  The  constant  part 
of  T  is  2  x  2,  real  skew-symmetric  and  can  thus 
be  synthesized  by  a  single  gyrator. 

This  -completes  the  proof  of  Theorem  19. 
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We  are  now  able  to  prove  the  following: 

Theorem  22  Let  Z  be  a  PR  matrix  of  degree  d .  Then 
the  synthesis  of  Z,  if  carried  out  by  the  hybrid 
matrix  method  described  in  Theorems  18,  19,  20 
and  21  will  result  in  a  lossless  buffer  containing 
d  reactors,  at  most  d  gyrators  and  a  termination 
of  resistors  and  gyrators,  plus  ideal  transformers. 

A 

Proof.  Assume  a  PR  matrix  Z  to  be  synthesized. 

The-  removal  of  an  entire  j-axis  pole  at  any  time, 
will  reduce  the  degree  by  precisely  the  degree  of 
the  pole  matrix,  which  appears  in  the  buffer. 

Suppose  a  Case  A  cycle  is  called  for.  Let 

Let  6Z  =  m.  Hence  6Z  =  m  +  1  since 
of  Theorem  18  is  1. 


0  -<  CD  <  oo 

o 


fb_j _ 4>_ 


<DJ 


Now  5Z"  =  2  and  5Z'  =  m  -  1.  The  absorp¬ 
tion  of  -b  into  the  perfectly  coupled  trans¬ 
former  requirsu  one  reactor  from  Z"  and  one  from 
Z' .  This  transformer  then  has  degree  1.  Thus 
the  buffer  has  degree  2,  one  due  to  the  trans¬ 
former  and  one  remaining  in  Z".  The  degree  of 
the  termination  is  m  -  2.  Thus  the  cycle  is  mini¬ 


mal. 
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Suppose  o)q  =  0  or  co.  Then  no  Brune  im¬ 
pedance  is  added.  5Z"  =  1  and  5Z*  =  m  -  1  and 
again  the  method  is  minimal. 

Suppose  a  Case  B  cycle  is  called  for.  Let 
6Z  =  m.  Then  6Z  =  5 (Z  +  B)  =  m  +  1,  since  6B  =  1 
and  B  and  Z  contain  no  common  poles . 

6Z"  =  2.  Hence  6Z'  =  m  -  1.  Both  reactors 
in  Z"  are  used  to  absorb  -B,  resulting  in  a  6 
terminal  reactive  buffer  of  degree  2.  Z'  con¬ 
tributed  a  reactive  matrix  of  degree  1  to  the 
buffer.  Thus  the  degree  of  the  termination  is 
m  -  2.  Thus  in  every  case  the  method  is  minimal. 

In  both  Cases  A  and  B,  Z"  contains  at  most 
1  gyrator.  In  case  B,  the  Brune  section  B(s) 
adds  at  most  1  gyrator  to  the  buffer.  Hence  in 
every  case,  the  degree  is  reduced  by  2  for  at 
most  2  gyrators  in  the  buffer. 

This  completes  the  proof  of  Theorem  22. 

In  the  following  pages  the  entire  synthesis 
procedure  is  illustrated  in  flow  chart  form.  The 
v  arious  steps  are  in  accord  with  the  lemmas  and 
theorems  of  this  chapter. 


Realize  with  R, G  and 
transformers 


Zjj  singular?  J) 


^  identically  sing 
ular  on  Re  s  =  0? 


Remove  resistance 


Re  Z  singular  at  jo>_? 


pply  congruence  trans-  Apply  congruence  trans¬ 
formation  to  give  formation  to  give 

Z  =  DtZ  D  and  (jcoo)  =  0  Z  =  D^Z  D  and  ) 

11  singular  11 


Cg 

identically  zero? 


Remove  gyrators 
and  realize  them 

^  is  the  remainder 


0  or  :o? 


Case  A 
Entry  2 


Case  A 
Entrv  1 


Case  B 
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Return  to  1 
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Example  of  Case  B  Synthesis 


In  order  to  demonstrate  the  Case  B  synthesis 
procedure,  we  have  selected  a  matrix  which  already 
has  the  property  that  the  hermitian  part  of  its 
first  principal  2x2  submatrix  is  singular  on  the 
j-axis.  Thus  the  initial  removal  of  resistance 
is  not  necessary  in  this  case.  Clearly  no  gener¬ 
ality  has  been  lost  by  selecting  such  a  pre¬ 
conditioned  matrix. 

We  wish  to  synthesize  Z  = 


37s 


(s 


^  +  43s  +  10 
+  1)  (s  +  2) 


6 


s  +  1 
s  +  2 


13s  + 
s  + 


6 


s  +  2 


S  +  1  ~  S  4-  1 

s  +  2  s  +  2 


11s  -  2 
s  +  2 


s  +  1  5s  +  6 

s  +  2  s  +  2 


Computing  Z  at  s  =  j  we  obtain 
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51  +  j  6  2  6(3  +  j)  41  +  j 12 


\ 


«<i> 


6(1  +  j2)  3  +  j  2(3  +  j) 


7  +  j  24  2(3  +  j) 


/  51  12- j3  24- j6^ 


1 

5 


12+j3  3 


^24+j6 


6  1+  | 


17 


62 


9+j6 


17  +  j4 / 
9-j6  18- j  17^ 

1  2 


^8+j 


17  2 


ZH^  +  ZSH^^  * 


Hence 


4  / 


z„  (j) 


H 


11 


1 

5 


51 


12+j3 


12-j3^ 


is  singular,  and  so  for  this  case  o>  =  1, 


Brume  Section 

In  order  to  make  Z^^j)  singular  it  is  neces¬ 
sary  to  add  a  Brune  section  which  cancels  (j)  . 

Ox*  ^ 
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In  Lemma  9  an  existence  theorem  is  given  which 

A t 

shows  that  a  matrix  B  exists  such  that  +  B 

is  singular  at  s  =  j.  In  what  follows  the  matrix 

B  is  obtained  in  a  slightly  different  fashion,  one 

which  is  perhaps  computationally  more  convenient. 

B(s)  will  be  so  chosen  that  B(j)  cancels  the 

matrix  Z.,„  (j).  To  gibe  B(s)  its  required  prop- 

SH11 

erties  we  also  add  an  undetermined  amount  of  the 

matrix  Z„  {')  .  The  determination  of  the  multir 
Hli  ' 

plier  finally  fixes  B  and  is  done  i  '  such  a  way 

i 

that  B  has  a  rank  1  pole  at  oo . 


j  62  9- j6  ^ 

f  51  12— j  3^ 

B(j)  =  -  \ 

-  j  ^ 

J  5 

^9+j6  X  1 

\12+j3  3  / 

where  X  is  chosen  so  that  B(s)  has  a  symmetric 
portion  comprised  of  a  rank  1  pole  at  oo,  plus 
a  constant  skew-symmetric  matrix.  The  symmetric 
part  of  B(j)  is 


B  ( j)  +  B(j), 


/62+51X  9+12x\ 


1 

5 


2 


9+12X 


1+3X 
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which  must 

be  made 

singular 

by 

the  choice 

of  X 

i.e.  (153X2+237X+62)  -  (144X2+216X+81)  = 

0 

This  gives 

X  =  : 

-7  ±  5/51 

6 

Using  X  =  - 

•7  -  5v/IT 
6 

gives 

/ 17  v5-l 
2 

1+2  Js\ 

< 

0 

v'S-A 

2 

B  ( j  )  =  j 

^1+2  5 

1+/5 

2  1 

+ 

(v/5-1) 

r  2 

7 

i.e.  we  add  to  Z  the  PR  matrix 

• 

(  17\/5-l 
2 

l+2\/5^ 

/  0 

/5-l\ 
2  * 

B  (s)  =  s 

^  1+2  \/5 

l+v^5  1 

2  / 

+ 

[_  h/s-1) 

\  2 

J 

Letting  Z  +  B  =  Z  we  obtain 

51 

12“  j 

Zn(j) 

■  1 1 

♦  i  1^)1 

^ 12+j  3 

3 

which  is  singular  as  required.  Z(s)  follows. 


143 


! 


I 

L 
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Z^  is  singular  at  s  =  j,  and  so  will 

have  a  pole  at  s  =  j .  Let 


u  -  h  \  zil  <  1 


“S12) 


Then  we  can  extract  the.  pole  matrix  from  U.  The 
residue  matrix  in  z”?-  at  s  =  j  is 


5)ZU 


3-75 

4 


Thus  the  residue  matrix  of  the  pole  in  U  at  s  =  j 


V  =  (s  -  j)U 


— ^  0  1 


1  -4+jj  1  0 .  - 

-4- j  17  0  1 


7-H24  6+i2 

15  5 


41+j 12 
5 

6+j  2 
5 
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3-y/5  -4- j 

4 


-4+j 


-3- j  2 


10+jll 


-3+j  2  10-jll 


which  is  of  rank  1. 

The  pole  at  s  =  -j  in  U  has  a  conjugate  residue  V* 
The  entire  pole  matrix  removable  from  *-s 


thus 


P(s)  =_l  +  41 

s-j  s+3 


- (8s+2) 


- (8s-2) 


- (6s+4) 


20s+22 


2 

s  +  1 


-  (6s-4 )' 


20s-22 


Letting  Z"  (s)  =  2(P)  we  obtain  z" 


/  34s 


1 

3^7? 


8s+2  11  (3-s/5 


8s-2 


2  (3~y/5) 


-ll(3-v/5)  -2(3-^'/ 
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l11 


4 

1 

0 


+ 


3±j5 

1 


0 


-2 


which  is  PR  of  degree  2, 
We  now  form  T 2^ 
following  page.  Q  =  ?  ^  2 
page  following  f ^  2(Z). 


with  1  gyrator. 
which  is  given  on  the 
(Z)  -  P  is  given  on  the 
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OKS 


-149- 


Then  Z' (s) 


ri<2(Q)  - 


fs  (62+29  V5) 
+  35+15/5 


i 

4 


(37+17  J5) 
2 

25+11  & 


I 

I 

I 

4 

i 


=  s 


(37+17/5)  1 

2  1 

17+7/5  1 

2  ! 

_  (H+5/5) 

8  2 

.  7+3 /5 

+  2 

4+5/5  1 

1+/5 

^62+29/5 

37+17/5 ' 
2 

37+17/5 

11+5/5 

2 

2 

\  0 

0 

/35+15/5 

21+9/5" 

4 

f 

21+9/5 

7+3/5" 

i  4 

2 

\3+/5 

0 
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2+^5 

2 


-2-Ts 


V 


2 

1+4  '/§’ 


-i-4  Til 
-l-v/5 


1+TiT 


The  first  of  these  three  terms  is  a  rank  1  pole 
at  co  and  can  be  realized  by  a  single  inductor 
and  ideal  transformers.  The  second  term  is 
purely  resistive  and  can  be  realized  by  3  resis¬ 
tors  and  ideal  transformers,  and  the  third  term 
can  be  realized  by  a  single  gyrator  and  ideal 
transformers. 

The  rank  1  pole  at  oo  in  Z*  is  called  M. 

It  is  caused  by  the  addition  of  B  to  Z, 

Prom  Z"(s)  we  obtain  a  series  pole  at  oo 
given  by 

f  17  4  0  ! 

4  10 

0  0  0 

The  connection  sequence  of  Z',  Z",  M,  N  and  -B 
is  shown  in  Figures  29,  30,  31  and  32. 


N 


s 


3+Ts 


The  matrices  -B,  N  and  M  of  Figure  32  are  com¬ 
bined  by  the  method  given  in  the  proof  of  Theorem 
21  into  a  matrix 


m  — 

j  -B  +  N  1 

I 

1 

N  ^ 

1  — 

l  ■  : 

N  + 

The 

reactive 

part  of  T 

is 

/  26 

5 

51+17/5 

2 

12+4/5  \ 
2  1 

f 

5 

1 

12+4/5 

2 

3+/5 

2 

s 

5 1+17/5 

12+4/5 

175+75/5 

49+21/5 

i  2 

2 

2 

2 

\  12+4/5 

3+J5 

49+21/5 

14+6/5  J 

\  2 

2 

2 

2  / 

which  is  rank  2  and  may  be  realized  by  con¬ 
gruence  transformer  and  2  inductors.  The  final 
realization  for  Z  is  shown  in  Figure  33. 


APPENDIX  I 


The  PR  Property 

A.  THE  IMPEDANCE  MATRIX. 

Definition  1.1  A  network  is  said  to  be 
passive  if  there  are  no  energy  sources  within 
that  network. 

Suppose  a  network  is  excited  by  an  external 
energy  source.  Any  increase  of  the  energy  with¬ 
in  the  netwoik  must  be  obtained  from  that  external 
source.  We  thus  have 

Axiom  1.1  A  passive  network  cannot  posess  an 
increasing  amount  of  energy  when  all  external 
forcing  functions  are  zero. 

Suppose  that  a  passive  network  is  excited 
by  an  external  forcing  function,  which  is  periodic 
(period  T) .  Once  a  steady  state  has  been  achieved 
the  transfer  of  average  energy  can  only  be  from 
the  exciting  source  into  the  network.  We  state 
this  as 

Axiom  1.2  The  average  energy  supplied  to  a  pas¬ 
sive  network  by  a  periodic  source  is,  in  the 
steady  state,  non-negative. 


-152- 


-153- 


Let  e (t)  be  the  instantaneous  energy  enter¬ 
ing  a  passive  network  and  p(t)  the  instantaneous 
power.  Then 


p<« 

Let  the  source  be  periodic  and  let  a  steady  state 
exist.  Then 


(n+1)  T 

P  «  ^  /  p  (t)  dt 

ave  T  J  ^ 

nT 

=  e((n+l)T)  -  e  (nT)  ]  . 

By  Axiom  1.2,  e((n+l)T)  -  e(nT).  Hence  we  have 
the  following 

Theorem  1 . 1  The  average  power  entering  a  passive 
network  from  a  periodic  source  is,  in  the  steady 
state,  non-negative. 

Suppose  that  there  are  n  ports  at  which 
forcing  functions  may  be  applied,  and  let  the 
responses  be  measured  at  those  n  ports.  Let  the 
network  consist  of  finite  numbers  of  resistors, 
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capacitors,  gyrators,  ideal  transformers,  self 
and  mutual  inductors,  all  of  which  are  linear. 
Let  the  forcing  functions  at  the  n  ports  be  cur¬ 
rents,  and  the  responses  be  voltages. 

Then  v^(t),  the  voltage  at  the  k-th  port 
is  related  to  ir(t),  the  current  entering  the 
r-th  port  by  a  linear  differential  equation  of 
the  form 


(1,1) 


P 

Z 

m=0 


%A) 


(t)  = 


q 

z 

v=0 


£  DVi  (t) 
rv  r 


where  D  means  ^  . 

The  coefficients  a.  and  B  are  obtained 

km  rv 

by  the  application  of  Kirchhoff's  laws  to  the 
network.  Since  a  network  comprised  of  a  finite 
number  of  elements  can  have  only  a  finite  number 
of  loops  and  nodes,  it  follows  that  p  and  q  are 
finite.  Moreover,  if  all  of  the  elements  in  the 
network  are  physically  realizable,  then  their 
values  are  given  by  real  numbers.  The  applica¬ 
tion  of  Kirchhoff's  equations  will  encounter  only 
ceal  numbers  and  so  all  of  the  and  will 

be  real. 


"T 
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Suppose  we  restrict  the  forcing  functions 
ir(t)  to  be  in  the  class  of  functions  which  are 
Laplace  transformable.  Assuming  no  stored  energy 
within  the  network  at  t  =  0,  we  obtain,  by  trans¬ 
forming  both  sides  of  Equation  1.1, 


where  s  is  the  Laplace  complex  variable. 
Define 


z.  is  a  rational  function  in  s  with  real  coef- 
kr 

ficients. 

By  the  assumption  of  linearity  .if  currents 
are  simultaneously  applied  to  all  n  ports, 


n 

V.  (s)  =  2  z^r (s)  ^ 

r=l 


k=l, 2. . .n 


This  gives  the  general  equation  for  the  network  as 


V (s)|  =  Z(s)I  (s)J  , 

where  V(s)[  is  the  vector  of  Laplace  transforms  of 
the  port  voltages  and  I(s){  the  vector  of  port 
currents.  Z(s)  is  an  n  x  n  matrix  called  the 
driving  point  impedance  matrix.  We  have  thus 
proved  the  following 

Theorem  1 . 2  Let  Z(s)  be  the  driving  point  im¬ 
pedance  matrix  of  a  passive  linear  finite  lumped 
(PLFL)  electrical  network.  Then  the  elements  of 
Z  ( s )  are  rational  functions  of  s  with  real  coef- 
ficients. 


B.  CONSTRAINTS  ON  THE  POLES  OF  z^y 

Let  all  currents  be  zero,  except  for  i  (t) . 
Suppose  that  ir(t)  is  a  rectangular  pulse  of 
finite  duration  starting  at  t  =  0,  and  such  that 
after  t  -  t^,  ir(t)  =  0.  If  zkr(s)  posesses 
poles  in  Re  s  >  0,  then  it  can  be  shown  that 
after  t  *=  t^,  when  all  forcing  functions  are 
zero,  v^(t)  increases  without  bound.  Since  this 
implies  an  increase  in  the  internal  energy  of 
the  system,  a  violation  of  Axiom  1.1  is 
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C.  ENERGY  CONSIDERATIONS 

Suppose  a  network  whose  impedance  matrix 
is  Z(s) ,  is  excited  by  a  current  vector  i(t){ 
where  ir(t)  is  entering  at  the  r-th  port.  Let 
i(t)\  be  of  the  form 
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where  and  are  all  real,  and  where 

o>o  /  0  or  oo .  Note  that  all  components  of 
i(t)|  are  sinusoidal  with  the  same  frequency. 


Phases  cx^  and  amplitudes  m^  are  completely 
arbitrary. 


i(t){  may  be  written 
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where  a\  and  bl  are  real  n-vectors.  Then 

i(t)l  =  Re{ejCUot  IU- 

Let  ( t )|  «  e*,CU°tll'.  Then  letting  H(s)l 

be  the  Laplace  transform  of  rj  ( t ){,  we  have 


H(s)| 


If  H(s)J  is  applied  to  the  ports  of  Z(s) 
the  corresponding  voltages  will  be 


N(s)}  =  Z(s)H(s)| 


I* 


Assume  that  ja)Q  is  not  a  pole  of  Z(s). 

By  Theorem  1.2,  j-axis  poles  in  z..(s)  are 

X  J 

simple,  and  poles  in  the  right  half  plane  do  not 
ekist.  Hence  letting  n(t)|  be  the  inverse 
Laplace  transform  of  N(s)|  we  obtain 
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n(tH  =  ztjaye^I*  +  Z  .  s-y  I| 


+  2  F„  (t)e  ^  +^ar^t  ij 


where  is  a  residue  matrix  at  a  possible 
j-axis  pole  of  Z(s)  and  where  Fr(t)  is  a  matrix 
of  polynomials  in  t  of  order  one  less  than  the 
multiplicity  of  the  r-th  pole  of  ^  (s)  in  the 
left  half  plane  at  s  =  +  jco^,  (o^c  * 

In  the  steady  state,  (t  very  large) 


n(t)(  -  Z(ja>o)e*B-t«  +  Z  ,  .  %)  ll 

n(t)|  has  a  real  and  an  imaginary  part.  The 
response  to  the  real  part  of  7)(t){  must  be  the 
real  part  of  n  (t)f  . 

Hence  if  i(t)f  =  Re  (e^CU°t  if)  we  obtain 

W  eja>ot 

v(tH  =  Re[Z(ja>o)eW  l|  +  2  ,  x|  ) 


The  instantaneous  steady  state  power  into  the 


network  is 
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n 

p  (t)  =  Z  i  (t)v  (t) 

*ss'  '  ,  mv  '  m'  ' 

m=l 

=  Re(e^CU°t'l)  Re(Z(  j^o)e^U)°t  l|  } 

W  e j<U-t 

+  Retell)  Re[Z  (* )  I|) 

k  J  '  k  o' 

Recall  that  1^  =  a|  4*  jbj .  Let 

Z(jo)o)I|  =  E{  =  c*  +  jdl 

and 

WkxV 

j  t<\  -  <bo)  "  ekt  +  jfkl- 

Then 

p  (t)  =  Re(  (a  +  Re{  (cl  +  jdl)e^°#t} 

O  b 

+  Re(  (a*  +  jb)e^ot)  Re(Z(ek*  +  jfkf  )e^ot  ) 

=  acl  cos  co  t  +  bdi  sin  o>  t 
o  o 

-  (adl  +  bcl  )sin  cn  t  cos  cn  t 

o  o 

+  (iT  cos  cn^t  -  "b  sin  )S(ev.}cos  o^t  -  fkJ  sin  a^t) 
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The  average  power  will  be 
0)  2ir/“o 

Pave  =  "2F  /  Pss(t>dt 

o 

=  -^(acl  +  bdl ) 


Now, 

\  Re  (l*Ej )  =  Re[(r  -  jb)  (c*  +  jdj  )] 

=  -^(acl  +  bdi )  . 

Thus 

Pave  =  1  Re  <r*E*  >  ' 

But  E|  =  Z(jcoo)iJ.  Hence 

Pave  =  4  Re<r*Z(jfflo)l|). 

But  by  Theorem  1.1,  Pave  ^  o.  Hence  we  must 
have,  for  not  a  pole  of  Z(s)  and 
for  0  -=•  co  -<  oo , 


Re  (l*Z  ( jo3Q)  i|)  ^  0. 


I  / 


vx 


•w 


>4* 


♦ 
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If  oj  =  0,  t)  (t){  =  ij  which  is  real  and 

n(t)|  =  Z(0)lj  which  is  also  real.  Hence 

pss(t)  =  Re (I)  Re ( Z ( 0 ) 1 1  ) 

+  terms  due  to  possible 

j-axis  poles  in  z.  .  (s) 

13 

=  Re  a  Re  cl 

=  acl 

Hence  P  =  "act  ,  and 
ave 

Re(I*EO  =  Re  (act)  =  act 
Thus  as  before,  for  0  not  a  pole  of  z^(s), 

Re (I*Z (0) l| )  ^  0. 

If  cdq  =  00  we  apply  the  transformation 
s'  =  -jjj-.  Then  s  =  00  maps  into  s'  =  0.  An 
inductor  sL  is  now  replaced  by  a  capacitor 

— ^ —  and  a  capacitor  <—7  is  replaced  by  an 
s '  (~)  SC 

inductor  s’ (”) .  Clearly  this  is  a  passive  net- 

work  and  as  before  P  -  0,  which  implies 

ave 
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Re  (I*Z (s)  I|  )  *  O  at  s'  =0 
We  have  thus  proved 

Theorem  1.4  If  Z(s)  is  the  driving  point  imped- 
ance  matrix  of  a  PLFL  network  and  joa  is  not  a 
pole  of  Z,  then  for  any  complex  vector  x, 

Re(x*Z(jo>o)x})  ^  0. 

Corollary  1.4.1  The  hermitian  part  of  Z(s)  is 
non-negative  definite  everywhere  on  the  j-axis 
where  Z  has  no  poles. 

Proof. 

Re(x*Z(s)xJ)  =  x*Z  (s)xj  for  all  3c 

H 

Corollary  1.4.2  If  Z(s)  is  the  impedance 

matrix  of  a  lossless  network  (i.e.  without  res- 

istors)  then  ZT,  (s)  is  null  everywhere  on  the 
H 

j-axis  where  Z  has  no  poles. 

Proof.  Since  the  average  power  delivered  to 
the  network  is  zero, 

Re  (x*Z(ju>Q)xJ  )  =  0 


for  all  x  . 
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D.  THE  ASSOCIATE  FUNCTIONS 

Definition  1.2  Let  £(s)  be  a  rational  function 
of  the  complex  variable  s,  with  real  coefficients. 
Then  f (s)  is  said  to  be  a  real  rational  function 
of  s . 

Definition  1.3  Let  F(s)  be  a  matrix  whose  ele¬ 
ments  are  rational  functions  of  s.  Then  F  is 
said  to  be  a  rational  matrix  function  of  s. 

F  may  be  a  rational  matrix  function  which  is 
any  or  all  of  square,  real,  symmetric  etc. 
Definition  1.4  Let  F  be  a  square  real  rational 
n  x  n  matrix  function  and  let  ~x  be  a  constant 
complex  n-vector.  Then  f(s)  =  x*Fxl  is  called 
an  associate  function  of  F.  Note  that  for  any 
F  there  are  arbitrarily  many  associate  functions. 

If  F  is  symmetric  then  f  is  a  real  rational 
function  for  any  3c  ,  If  3T  is  real  then  f  is 
a  real  rational  function  even  if  F  is  not  sym¬ 
metric. 

In  general  however,  f  is  a  rational  function 
with  complex  coefficients. 

Theorem  1.5  Let  Z(s)  be  the  matrix  of  a  PLFL 
electrical  network,  Let  z(s)  be  an  associate 


iiaig.ti.Qn-Qf  z.  Haai 


i )  Re  z  ( s )  -  0  for  Re  s  =  0 
i i )  z ( s )  has  no  poles  in  Re  s  >  0 
iii)  I_f  z(s)  has  poles  on  Re  s  =  0 

then  they  are  of  single  multiplicity, 
and  have  positive  real  residues. 


Proof. 

Proposition  i) 
Theorem  1.4. 

By  Theorem  1.3, 
in  Re  s  :>  0.  Now  z 


follows  immediately  from 


none  of  the  z. .  have  poles 
n  n  x3 

=  2  2  z  .  .x.  *x  .  ,  and 

i=l  j=l  1  3 


so  z(s)  cannot  have  poles  in  Re  s  >  0.  Hence 
ii)  is  proved. 

By  Theorem  1.3,  if  any  of  the  z ^  have  poles 
on  Re  s  =  0,  then  they  are  of  simple  multiplicity. 
Hence  z(s)  can  have  poles  on  Re  s  =  0  only  if 
they  are  simple. 

Suppose  z(s)  has  a  pole  at  Let 

j0 

s-jcdo  =  peJ  . 


Then  for  very  small  p, 
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IQ 

where  Ke  is  the  residue  at  the  pole.  Since 
Re  z(s)  ^  0  on  s  =  jio  (by  i)  we  have 

cos  (a  +  9)  s?  0  for  6  =  7r/2,  -tt/2. 


But  this  can  only  be  true  if  a  =  0,  i.e.  z(s) 
has  a  positive  real  residue  at  a  j-axis  pole. 
This  completes  iii)  and  proves  the  theorem. 
Corollary  1.5.1  If  Z,  the  matrix  of  a  PLFL 
network,  has  a  pole  on  the  j-axis  then  the 
residue  is  a  non-negative  definite  hermitian 
matrix. 


Proof. 


Let  z(s)  -  x*Z(s)x|  .  If  Z(s)  has  a  pole  at 

be  the 


s  =  j0)Q,  let  W  =  (s  -  ja)Q)Z(s) 


residue  matrix,  and  let 


s=joo 


w  =  (s  -  jo>o)z(s) 


be  the 


s=jai 


residue  in  z(s).  Then 


w  -  x*Wxj  . 


But  if  w  if  non-zero,  then  by  iii)  of  Theorem 
1.5,  it  must  be  positive  and  real.  Hence 
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x*Wx|  -  0. 

This  m*  ans  that  W„  is  ncn-negative  definite  and 

Wotl  is  null.  Hence  W  is  non-negative  definite 
on 

hermitian,  which  proves  the  corrolary. 


E.  FURTHER  IMPLICATIONS 

Lemma  1.1  Excluding  possible  j-axis  poles,  let 
f(s) 

i)  be  requalar  on  Re  s  =  0 

ii)  have  no  poles  in  Re  s  >  0 

iii)  be  such  that  Re  f(s)  -  0  on 
Re  s  =  0. 

Then  Re  f ( s )  -  0  for  Re  s  ^  0. 

Proof. 

Let  C  be  the  closed  curve  consisting  of  the 
j-axis  except  at  poles  of  f(s)  where  C  is  a 
small  semi-circle  of  radius  p  in  the  right  half 
plane.  Let  R  be  the  region  enclosed  by  C,  in 
the  right  half  plane. 

By  the  Principle  of  the  Minimum,  (well 
known  in  functions  of  a  complex  variable) 


0  -  Re  f  { 8 ) 


*  Re  !:  (s) 


s  e  c 


s  £  R 
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Since  p  can  be  taken  arbitrarily  small  it  follows 
that,  except  at  possible  j-axis  poles, 

Re  f(s)  ^  0  for  Re  s  ^  0. 

This  proves  the  lemma. 

Corollary  1.1.1  If  z(s)  is  an  associate  function 
of  Z(s),  the  matrix  of  a  PLFL  network,  then 

Re  z(s)  -  0  for  Re  s  -  0. 


Proof. 

By  Theorem  1.6,  all  conditions  of  Lemma  1.1 
are  satisfied  for  z(s). 

We  can  now  prove  the  following 
Theorem  1 . 6  Let  Z ( s )  be  the  matrix  of  a  PLFL 
network.  Then  ZT,  (s)  is  either  hermitian  positive 

,  -  -  J-J 

definite  in  Re  s  >  0  or  it  is  identically  null 
for  all  s. 

Proof. 

Let  z(s)  be  an  associate  function  of  Z(s). 

Suppose  Re  z(s)  =  0  at  an  isolated  point  s  =  s^ 
in  the  right  half  plane,  and  suppose  that  z(s)' 
is  not  identically  zero.  i.e. 

z(so)  =  j3  (3  real) 


Then 


g(s) 


where  g(s)  has  poles  at  sq  of  order  n-1  or  less. 
_  j  0 

Let  s  -  sq  =  pe  .  Then  for  p  sufficiently 
small, 


^'zTs)1-  j$>  *  f  cos(a  -  n®> 


Now 


Re  ( 


z 


Re  z  ( s ) 

(Re  z  ( s ) )  2  +  [im  z(s)  -  p)  2 


But  by  the  corollary  to  Lemma  1.1,  Re  z(s)  -  0 

in  Re  s  ^  0.  in  particular,  this  is  true  on 
j  9 

s  -  sq  =  pe J  .  Hence 

cos  (a  -  n0)  5:  0  for  0^6*  2w . 


Clearly  then  n  =  0.  But  in  that  case  Re  z(s) 
does  not  have  a  zero  at  an  isolated  point  in 
Re  s  >•  0,  which  contradicts  the  assumption  that 
it  does. 

Since  Re  z(s)  ^  0  in  Re  s  ^  0,  either 
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Re  z(s)  >*  0  in  Re  s  >  0  or  Re  z(s)  =  0  for  all  s. 

Since  this  is  true  for  any  associate  function  of 

Z(s)  it  follows  that  either  Z  is  hermitian 

H 

positive  definite  in  Re  s  >  0  or  it  is  null  for 
all  s.  This  completes  the  proof  of  the  theorem. 

P.  THE  PR  PROPERTY 

Definition  1.5  A  square  matrix  function  F(s) 
is  said  to  be  positive  real  (PR)  if  I  and  II 
are  satisfied. 

I  The  matrix  elements  f. .  are  rational 
in  s  with  real  coefficients 
II  For  any  complex  vector  jT  , 

Re  (x*Fx|  )  -  0  for  Re  s  ^  0. 

Theorem  1.7  Let  F(s)  be  a  square  real  rational 
matrix  function.  Let  f ( s )  be  an  associate  func¬ 
tion.  Then  F  ( s )  _is_  PR  if  and  only  if 

lla.  Re  f(s)  -  0  for  Re  s  =  0 

l lb .  f ( s )  has  no  poles  in  Re  s  >  0 
lie.  For  Re  s  =  0,  poles  of  f(s)  are 

simple  and  have  non-negative  residues. 

Proof. 

Suppose  F(s)  is  PR.  Then  for  any  vector  x  , 
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Re  f (s)  =  Re(x*Fx|  )  -  0,  on  Re 


Suppose  f(s)  has  a  pole  in  Re  s  >  0. 


Ke 


Da 


f  (s)  = 


<s  -  So> 


n 


+  g  (s) , 


s  =  0 

i .  e . 


where  g(s)  has  poles  of  order  n-1  or  less  at  sq. 

•  Q 

Let  s  -  sq  =  pe-1  .  Then  for  very  small  p, 


f(s) 


K  j  (a  -  n0) 

P 


and  Re  f(s)  =  §  cos  (a  -  n0) ,  (o  -  0  -  2ir) . 

r 

But  this  is  non-negative  which  can  only  be  true 
if  n  =  0.  Hence  f(s)  cannot  have  poles  in 
Re  s  >*  0. 

Suppose  f(s)  has  a  pole  on  the  j-axis. 
Then,  following  the  proof  of  proposition  iii) 
of  Theorem  1.5,  such  poles  of  f(s)  are  simple 
and  have  non-negative  residues. 

This  completes  the  proof  of  the  theorem. 
Corollary  1.7.1  Let  F ( s )  be  a  PR  matrix.  If 
F(s)  has  a  j-axis  pole  then  the  residue  matrix 
is  hermitian  non-negative  definite. 
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Theorem  1.8  Let  Z(s)  be  the  impedance  matrix 
of  a  PLFL  network.  Then  Z(s)  _is  PR. 

Proof. 

By  Theorem  1.2  the  z. .  are  rational  real  functions 

of  3. 

By  Theorem  1.4,  for  any  vector  x  , 

Re  (3T*Zx|  )  ^  0  on  Re  s  =  0. 


By  Theorem  1.6,  either 

Re  (x*Zx|  )  >  0 
or  it  ,is  identically  zero. 

Re  (x*ZxJ  )  0 


in  Re  s  >  0 

Hence 

in  Re  s  ^  0. 


This  proves  the  theorem  and  completes  the  appendix. 


APPENDIX  II 


Diagonalization  of  Skew-symmetric  Matrices. 

Let  A  be  a  real  skew- symmetric  matrix  of 
order  2n.  Then  there  exists  an  orthogonal 
matrix  V  such  that 


where  the  are  real  and  possibly  zero.  This 
is  proved  as  follows. 

Lemma  2 . 1  The  eigenvalues  of  a  real  skew  sym¬ 
metric  matrix  are  pure  imaginary,  or  zero. 
Proof. 

Let  A  be  real  skew  symmetric,  let  X  be  an 
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eigenvalue  with  associated  eigenvector  x. 
Then 


Ax{ 


Ax} . 1 


and 


Ax*t  =  A*x*J  ...  2. 

Hence  x*Axl  =  Ax*x{  (by  1)  and  X*AxJ  =  -A*x*xl 
by  2.  Since  x  is  not  the  null  vector 

A  =  -A*.  QED 

Lemma  2.2  Let  A  be  a  real  skew  symmetric  matrix, 
let  ip.  be  an  eigenvalue  (p.  real,  nonzero)  and  let 
x\  +  iyl  be  the  associated  eigenvector,  where  the 
vectors  xj  and  y|  are  real.  Then  x|  and  yj  are 
orthogonal. 

Proof. 


A  (xj  +  iyl  )  =  ip.  (xl  +  iyj  ) 


Thus 


Ax  l  =  -p.yl 


1 


and 
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AyJ  =  |ix| . 2. 

Hence  by  2,  y(AyH  =  -nxyl  ,  giving  uyxl  =  -u5ty*  . 
But  M- /  xj  and  yj  are  real,  nonzero. 

Thus  xyj  =  0.  QED 

Lemma  2.3  Let  A  be  a  real  skew-symmetric  n  x  n 
matrix,  i|x  a  nonzero  eigenvalue  with  associated 
eigenvector  x}  +  iyj  .  Let  M  be  a  real  ortho¬ 
gonal  matrix  whose  first  two  columns  are  x\  and  y| 
Then 

0  4  | 

I 
i 

4  0  | 

- 1  i 

*  !  An-2  / 


where  ^  is  a  real  skew  symmetric  matrix  of 
order  n-2. 

Proof. 


Let  M  =  (xj  y\ 

•  • 

.ml) 

nf 

be  an  orthogonal 

matrix.  Then 

AM  =  {  Axj 

Ayl 

Am3( 

.  .  .  Am  i  ) 
n*  ' 

=  (-4YI 

4*1 

Am3) 

.  .  .  Am  1 ) 

n* 

/ 

M  = 

\ 
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Now 


xAm^J  =  -rn^Ax^ 
= 

=  0. 


Hence 


Since  A  is  skew-symmetric,  so  is  V^A  V.  Hence 

so  is  A  _ 
n-2. 


QED 
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Theorem  2 . i  Let  A  be  a  real  skew-symmetric  matrix 
of  order  2n.  Then  there  exists  an  orthogonal 
matrix  V  such  that 


where  the  jju  are  real  and  some  possibly  zero. 
Proof. 

If  A  has  only  zero  eigenvalues  then  it  is  the 
null  matrix  and  we  take  V  as  the  identity  matrix, 
and  the  proof  is  complete.  If  A  has  a  nonzero 
eigenvalue,  let  it  be  called  ip^.  Then  by  Lemma 
2.3,  there  exists  an  orthogonal  matrix  M-^  such 
that 


T- 
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Suppose,  for  a  proof  by  induction,  that  there 

exists  an  orthoaonal  matrix  M  such  that 

r 


(We  call  this  the  r-fold  skew-diagortalization  of 

A).  If  A0/  v  is  the  null  matrix  the  rest  of 
2 (n-r) 

the  proof  is  trivial. 


% 
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let 


Dr+1A2 (n-r)Dr+l 


/ 


“M- 


r+1 


^r+1  ! 


0  i 


\ 


0 


\ 


0: 


A2 (n-r-l) j 


The  existence  of  is  guaranteed  by  Lenuna  2.3 


Hence  K  M  K  , 

r  +  1  r  r  r+-l 


But  M  Kr+^  is  a  product  of  orthogonal  matrices 
and  so  is  again  orthogonal.  Let  M  ^  =  M  K 
Then  we  have  shown  that  if  there  exists  an  ortho 
gonal  matrix  <>:.  which  is  supposed  to  skew-dia- 
gonalize  A  to  an  r-fold  level,  there  also  exists 
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an  orthogonal  matrix  Mr+^  which  skew-diagonalizes 
A  to  an  (r+l)-fold  level. 

But  exists.  Hence  Mr  exists  for  all  r 
positive  and  integral,  where  the  matrix  A  has 
the  nonzero  eigenvalues  ip^.  -  .  .  ipr  (r  -  n)  . 
Let  V  =  M^M2*  .  ,Mn.  Then  V  is  orthogonal  and 


where  some  of  the  \i ^  are  zero.  This  proves  the 
theorem. 

Corollary.  Let  A  be  real  skew-symmetric  of  Drder 
2n+l.  Then  there  exists  an  orthogonal  matrix  V 


such  that 
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Figure  5  -  Realization  of  Z(s)  of  Example  2 


gure 
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Figure  11  is  deleted. 


TRANSFORMER 
MATRIX:  D"1 

FIGURE  12”  CONDITIONING  Z  FOR  SYNTHESIS 


FIGURE  I3~  SYNTHESIS  OF  THE  FIRST  PART  OF  Z" 


FIGURE  14-SYNTHESIS  OF  SECOND  PART  OF  Z 
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PERFEGTLY  COUPLED 
TRANSFORMER 

TURNS  RATIO  MATRIX: 
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FIGURE  17-  SECOND  TERM  OF  Z"  (ALTERNATE  REALIZATION) 
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FIGURE  18-  ALTERNATE  REALIZ 


CONGRUENCE 

TRANSFORMER 


•4 


207 


FIGURE  28-  DECOMPOSITION  OF 
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FIGURE  27-  COMPENSATING  FOR  THE  BRUNE  SECTION 
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FIGURE  32-  Z  { S )  —  B(S) 
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FIGURE  34-  KYPERCUBE  OF  GYRATION  MATRICES 
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FIGURE  36—  HYPERCUBE  FOR  Z=S 
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